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Abstract 

A model is proposed to demonstrate that classical general relativity can emerge 
from loop quantum gravity, in a relational description of gravitational field in terms of 
coordinates given by matter. Local Dirac observables and coherent states are defined 
to explore physical content of the model. Expectation values of commutators between 
the observables for the coherent states recover the four-dimensional diffeomorphism 
algebra and the large-scale dynamics of the gravitational field relative to the matter 
coordinates. Both results conform with general relativity up to calculable corrections 
near singularities. 
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1 The upper figure pictures a generic part of an embedded graph, which is an 
embedded edge e.\ connecting the two nodes v<i and V3. The lower figures 
demonstrate the corresponding definitions of e, + , e ir , E31 and €21- In this 

V P ' 

case, e 2) i is the same as ei [J3] 

2 The left figure pictures an embedded graph 7, with one of its nodes labeled 
by n = 1 and the four edges connected to this node labeled by (1, j). The 
figures in (a) demonstrate the corresponding definitions of v%, ■ and e' ni p 
while the figures in (b) demonstrate the corresponding definitions of ej n ^ and 

ej [22] 

3 The lattice rectangular prism 7| C M 3 with the vertices represented by the 
black dots at the locations and the links represented by the lines be- 
tween two adjacent vertices. The lattice torus 7t orus for A can be obtained by 
identifying the opposite boundary faces of if. As an illustration for chapter 8, 
the construction also assigns the location X n of a vertex to the corresponding 
vertex v* of the lattice torus as its coordinate value [33] 
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In our dreams (writes Coleridge) images represent the sensations we think they 
cause ; we do not feel horror because we are threatened by a sphinx ; we dream a 

sphinx in order to explain the horror we feel. If this is so, how could a mere 
chronicle of its forms transmit the stupor, the exaltation, the alarm, the menace 
and the jubilance which made up the fabric of that dream that night ? 

- J. L. Borges, Ragnarok 



vii 



1 



1 Introduction 

Loop quantum gravity [4J [5J [6 J is a candidate quantum theory of gravity. The theory strives 
to construct the deep quantum structure of spacetime, while holding on to background inde- 
pendence that is the heart of general relativity. Following the approach of Dirac quantization 
in the Arnowitt-Deser-Misner (ADM) formalism [4][5J[6J, the theory's kinematical Hilbert 
space, called knot space, is rigorously constructed and is shown to be almost unique. Knot 
space describes spatial geometry in terms of quanta of areas and volumes, realizing the idea 
that the geometry of space should be discrete with Planck-sized units. 

Because of its novelty, the theory also faces new challenges [4] [5] [6] : based on the Planck- 
scale discrete structure of the knot states, the recovery of the observed smooth geometry has 
been a major issue for the theory; further, the background independence demands one to 
solve the Hamiltonian constraint upon knot space to obtain a physical Hilbert space, and 
also to construct diffeomorphism-invariant local Dirac observables. The difficulty in these 
tasks makes the semi-classical limit of the theory hard to obtain, and thus it is unknown 
whether loop quantum gravity describes the gravitational field at observed scales. However, 
the tasks are completed in a symmetrically reduced setting in loop quantum cosmology 
[31j[30j[28j[29j, minisuperspace models that carry the main features of loop quantum gravity. 
Moreover, it has been shown that loop quantum cosmology has correct semi-classical limits, 
with meaningful corrections in the regions where the classical theory breaks down [3Tj [30J . 
Therefore, loop quantum cosmology provides a valuable guidance for the full theory. 

Following this guidance, the model proposed in this thesis explores the semi-classical 
limit of loop quantum gravity. The model shows that the above tasks can be accomplished 
with a few explicitly specified assumptions. Further, the obtained semi-classical limit is 
also shown to reproduce general relativity. The system considered by the model contains 
both gravitational and matter fields, and both sectors are treated with the standard loop 
quantization procedure to obtain the kinematical Hilbert space for the system. To go beyond 
the kinematical level, three explicit assumptions are made: 1) that a modified Hamiltonian 
constraint operator is valid; 2) that group averaging method solves the modified Hamiltonian 
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constraint and gives the physical Hilbert space for the model; 3) that the matter back- 
reactions on the gravitational dynamics can be ignored in this context. 

The matter fields, with their negligible back reactions, provide internal spacetime coor- 
dinates for the local observables in the model. Using matter coordinates to describe local 
dynamics is a long existing idea in general relativity [19| [18J [17] , the implementation of which 
in the quantum theory leads to the local quantum observables [15j[16j[18j[17j. The inter- 
nal coordinates are provided by composite matter fields, and are also quantized. Since the 
coordinate fields are part of the fully coupled system, the observables are thus defined in a 
background independent manner and are diffeomorphism invariant. 

Subsequently, the appropriate coherent states in the model are defined to minimize the 
uncertainty of the gravitational local observables. The large scale physics is then obtained by 
explicit calculations based on the observables and the coherent states. First of all, the matter 
coordinates provide a natural correspondence between the expectation values of the observ- 
ables and a classical emergent gravitational fields. In terms of the emergent gravitational 
field, the quantum constraint algebra of the model reproduces the classical diffeomorphism 
algebra up to quantum corrections. Further, the emergent gravitational field satisfies clas- 
sical constraint equations up to quantum corrections. Moreover, the clock time derivatives 
of the emergent fields are shown to be given by the temporal translations generated by the 
quantum algebra. Finally, the equations of motion recover the classical gravitational dy- 
namics of general relativity up to quantum corrections. All of the correction terms for the 
semi-classical limits are small in large scales and nonsingular regions of spacetime. Therefore, 
the model's semi-classical limit reproduces general relativity, up to significant and calculable 
corrections where the classical theory breaks down. 

2 Loop Quantum Gravity 

In the first part of the thesis, I briefly introduce the basics of the loop quantum gravity, and 
describe how it couples to the matter fields. The purpose is mainly to explain the kinematics 
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of the theory, which captures quantum geometry of space background independently, and 
provides a foundation for the dynamics. 

2.1 Classical Hamiltonian Formalism 

Loop quantum gravity is the result of quantizing general relativity in a special Hamiltonian 
formulation, which uses triad and connection fields to describe gravity [lj [2 J [3J . 

The most common form of Hamiltonian general relativity uses metric and extrinsic 
curvature of space as phase space variables PQ[2J[3J. Instead of using the spatial metric 
<7a&(x), we can introduce triad fields {ef(x)}(x = (x,y,z);a = x,y,z;i = 1,2,3) consisting 
of three vector fields that are orthonormal according to the metric. In this correspon- 
dence, each set of triad fields specifies a spatial metric by q a b(x) = <5jj e a( x ) e b( x ) uniquely 
(e^(x)e"(x) = 5j). On the other hand, setting U?j(x) to be a local SO (3) transformation, 
e"(x) and ef(x) = C//(x)e"(x) would determine the same g &(x). Thus the use of the triad 
fields in place of the metric introduces an additional local SO (3) symmetry. 

The triad fields also serve as a spatial frame in space, such that we can write every spatial 
tensor using the triad fields as the basis. For instance, a vector field will be written as V l (x). 
In this frame, the spatial Levi-Civita connection will be written in the form: 



Moreover, the canonical conjugate momenta of e"(x) is the spatial extrinsic curvature field 
fQ(x), whose explicit form is given by: 



where C n denotes the Lie derivative with respect to the unit normal field n perpendicular 
to the spatial slices (according to the spacetime metric). There's no surprise that general 
relativity can be written in terms of the triad fields serving as an alternative basis. However, 




#j(x) = 6^K ab e%x) ; K ab (x) = C n q ab (x) 
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this new point of view leads to the discovery of Ashtekar formalism [2J [3 J , which reveals the 
striking similarities between general relativity and Yang-Mills gauge theory. 

The fundamental variables in Ashtekar formalism are related to the triad and spatial 
extrinsic curvature fields, through a canonical transformation 

E«(x) ee det(eK(x) 

(2.1) 

4(x)=rl(x) + 7 ^(x) 

where the real number 7 is called Immirzi parameter, and classically different values for 
the parameter correspond to different canonical transformations in (2.1). By construction, 
the fields Ef(x) are densitized triad fields, and the field A l a (x) are SO (3) gauge fields. The 
variables also satisfy the simple Poisson algebra, with the only non-vanishing bracket 



{4(x),£?(y)} = 87r(G/c 3 ) 7 ^(x,y) (2.2) 

where G is the Newton's constant. Using (2.2), the Hamiltonian general relativity can be 
reformulated in terms of the new canonnical conjugate variables A l a (x) and Ef(x). Details 
of the reformulation can be found in chapter 2 of [B] and chapter 3 of [12J. The resulting 
theory has the Hamiltonian 



R g (N, A, V) = H g (N) + G g (A) + M g {V) (2.3) 

where each of the three terms comes from smearing a constraint functional with an arbitrary 
Lagrangian multiplier (with the bars to emphasize the non-dynamical nature) over the spatial 
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manifold M 



H g (N)= [ d 3 xN(x)H g (x) = [ d *xN(x)^L 
Jm Jm v det E 



^ k F k ab + 2{l-^)K\ a Ki ] 



G g (A) = [ d 3 xA\x)G g ,(x) = [ d 3 xA l (d a E? + e tl k AiE a k )(x) 
Jm Jm 



(2.4) 



M g (V) = / d 3 xV a (x)M g<a (x) = 

Jm jm 



d 3 xV a (x) ( E\Fi b - l —t K i a G g , i ) (x) 



where K l a is now considered a function of the canonical variables A and E, as determined 
from (2.2) and (2.1). It is known that general relativity is a theory of pure constraints, and 
in this formalism the solutions to the theory are given by the constraint equations 



H g (x) = M gA (x) = G gji (x) = (2.5) 

The smeared Hamiltonian constraint H g (N), with any lapse function N(x), generates a dif- 
feomorphism on the fields that satisfy the on-shell condition (2.5) in the directions perpen- 
dicular to spatial slices. The smeared momentum constraint M g (V), with any shift function 
V a (x), generates a spatial diffeomorphism. Together, the smeared Hamiltonian and mo- 
mentum constraints generate spacetime diffeomorphisms. Additionally, the smeared Gauss 
constraint G g (A), with any A 1 , generates a local rotation of the triads. Clearly, the presence 
of Gauss constraint is due to the additional local SO (3) symmetry introduced along with 
the triad fields [5] [I] [6]. 

From here on, we replace the SO (3) symmetry group with SU(2) group, exploiting the 
fact that SU(2) group is the universal covering of SO (3) group and thus has the same Lie 
algebra. The constraints (2.4) and Poisson bracket (2.2) remain exactly the same under the 
replacement, so the physical content of the theory is unchanged. Under such illumination, 
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the similarities between general relativity and 577(2) Yang-Mills theory are revealed. The 
two theories share the same preliminary phase space parameterized by the fields Ef(x) and 
A l a (x), which represent (generalized) electric fields and magnetic potential fields in SU(2) 
Yang-Mills theory. Further, they are both SU (2) symmetric and they share exactly the same 
Gauss constraint. These similarities inspire the introduction of the loop representation from 
gauge theories into quantum gravity, and result in significant parallelism of the two in their 
quantization procedures. However, there is a major distinction between the theories that 
causes drastic departures of their physical Hilbert spaces. 

The difference lies in the background independence of general relativity. Treating space- 
time as a dynamical entity instead of a fixed physical background, general relativity uniquely 
possesses diffeomorphism symmetry. As a result, in contrast to the case of SU(2) Yang-Mills 
theory, the value of an SU (2) invariant field at a spacetime coordinate point (x, t) is not an 
observable, since it is not invariant under any diffeomorphism that moves the coordinate 
point. Further, the Hamiltonian of the theory is composed of the constraints. That means 
the evolution along a Hamiltonian flow is merely a gauge transformation, and any observable 
must be constant along the flow. To incorporate these features, the central theme of the 
quantization procedure toward loop quantum gravity is about achieving a quantum theory 
that is devoid of any fixed physical background. 

As a theory of pure symmetry, general relativity is governed by the algebra of diffeo- 
morphism and local SU (2) group. Denote [A, A'] to be the SU (2) commutator, and set 
[V,V'] a = V b d b V' a - V' b d b V a . In terms of the generators H g (N), M g (V) and G g (A) the 
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algebra is 



{Gg(A),G g (A')} = 8 7 r(G/c 3 ) 7 G 9 ([A,A / ]); {G g (A), M g {V)} = 8n(G/c 3 )jG g {CvA) 



{M g (V),M g (V')} = 87r(G/c 3 ) 7 M 9 ([y,n) 



{G ff (A), if ff (JV)} = 0; {M g (V), H g (N)} = 8tt(G /c^H^CyN) 



{H g (N),H g (N')} = 87t(G/c 3 ) 7 (M a (S) + G g (S a A a )) + ^ (r ?^ lG g ( [t "°' ,X ' E ''°" S '' 



^(G/c 3 



detSl 



(2.6) 

where £y denotes a Lie derivative and 

S a = (Nd b N' - N'd b N) T -i— : - (2.7) 

| det £/ 1 

Note that instead of structure constants, the algebra has structure functionals depending on 
the phase space variables. 



2.2 Quantization in Loop Representation 

Loop quantum gravity results from Dirac quantization of general relativity, in the Ashtekar 
formalism based on loop variables [5j[4j[6j. Explicitly background independent, the theory 
is non-perturbative, in contrast to the usual forms of particle theories. Originating in non- 
perturbative quantization of Yang-Mills theory, the loop representation is introduced to 
quantize gravitational fields in the theory. As a kinematical quantum space, knot space is 
rigorously constructed to solve Gauss and momentum constraints. The knot states in the 
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space describe spatial quantum geometry. The construction of the physical Hilbert space is 
an ongoing project involving quantizing and solving Hamiltonian constraint in (dual) knot 
space. 

2.2.1 Cylindrical Functions 

Based on the Ashtekar formalism given in the previous chapter, loop quantum gravity uses 
the connection fields A l a (x) as configuration variables to construct quantum states. To employ 
the loop representation we introduce a special basis, called cylindrical functions, for general 
wave functionals of the connection fields [5J [4 J [6J . 

Cylindrical functions are wave functionals that depend on the connection fields specifically 
through holonomies. A holonomy is a group-valued functional of connection fields defined 
with an oriented path e G M in the spatial manifold M. The bar indicates that e is 
embedded in M. Here, each holonomy will be SU(2) valued, and has an explicit spin j 
matrix representation as 



where k denotes SU(2) indices in the spin j representation, and V denotes path ordering. 
Under the transformations of the connection fields, a holonomy behaves simply under local 



where e(o) and e(i) denote the beginning and end point of e. Given a specific connection fields 

A l b (x), a holonomy specifies the parallel transportation of a spinor across the path e through 

the linear map (2.8) . Moreover, if we consistently assign group values to all holonomies 

along all paths in M, we uniquely specify a distributional connection fieldj^] Relying on this 

■'■Every smooth connection field on M gives a consistent assignment to the values for all holonomies. 
However, all possible consistent assignments of holonomy values, in general, correspond to non-smooth 
connection fields that are distributional. 




(2.8) 





(2.9) 



9 



fact, loop representation uses holonomies to build cylindrical functions for the basis. 

Define a graph 7 in M to be a set {e^} containing N e oriented paths, called edges, meeting 
at most at their end points. Again, the bars indicate that 7 is embedded in M. A cylindrical 
function t/j^j is specified by a graph 7 and a complex function / : SU (2) Ne — > C. Also, each 
cylindrical function is a wave functional of the connection fields through the relation 



The vector space spanned by cylindrical functions is denoted by Cyl, and carries an (over) 
complete representation of the wave functions. The space is endowed with inner products 
given by 



where dh is the normalized Haar measure of SU(2). The space Cyl with its inner product 
then serves as the primary space for the quantization to proceed. 

Recall that general relativity is given by the three constraints (2.4). Following the Dirac 
quantization approach, loop quantum gravity imposes the constraint equations (2.5) by 
demanding the physical states to be annihilated by the quantized constraints. We denote S 
the space obtained from imposing Gauss constraint on Cyl, and K the space obtained from 
imposing momentum constraint on S. The currently unfinished step is the construction of 
the physical Hilbert space H by imposing Hamiltonian constraint on K. In the following we 
will briefly describe the well-established part of the procedure that successfully leads to S 
and K. 



^j[A] = /(MeO, h(e 2 ), h(e 3 ), h(e Ne ))[A] 



(2.10) 




I II dh^nKex), He*), He 3 ), h(e N .))f(hft), h(e' 2 ), h(e' 3 ), h(e' N ,)) 



(2.11) 



eG7U7' 
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2.2.2 Space S of Spin Network States 

The first constraint to be imposed is the Gauss constraint, which generates local SU(2) 
transformations in the classical theory. The action of an arbitrary local SU(2) transformation 
U(x) on Cyl is directly given by Uu 

(2.12) 

MMei), h(e 2 ), h(e 3 ), h(e N J) = f{Uh{e x ), Uh{e 2 ), Uh(e 3 ), Uh(e Ne )) 

The quantized Gauss constraints are the generators of the finite transformations (2.12), thus 
the states satisfying the constraints should be invariant under the transformations. 

It turns out that the SU(2) invariant subspace of Cyl is spanned by the cylindrical 
functions called spin network states [5]|4][6j. Each spin network state is characterized by 
the following information (fig.l): 1) a closed embedded graph 7 in M, with its N e oriented 
edges {ej} connecting to its N v nodes {v n } without any free end; 2) a non-trivial SU(2) 
spin representation ji assigned to each edge; 3) a generalized Clebsch-Gordan coefficient 
(intertwiner) i n assigned to each node, with the exactly right indices to contract invariantly 
with the matrix indices of the holonomies along the edges meeting at v n . The explicit form 
of a spin network state S^^ u i n as a cylindrical function is 

{N v N e ^ 

(gK(g)/^%>) (2.13) 
n i ) 

where Inv{...} denotes the local SU(2) invariant tensorial product. Thus, (2.13) explicitly 
shows the invariance of Syj^ under the transformations (2.12). Note that the closure of 7 
is essential for the application of Inv{...}. With the inner products given by (2.11), the set 
of all spin network states {S^j it i n } provides a complete (but not over complete) orthonormal 
basis for S which solves Gauss constraint [S] [4| [6 J . 

Note that the orthonormal basis {S^j it i n } contains spin network states that differ only 
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by arbitrarily small changes of the embedded graphs. That means, in contrast to the case of 
Fock space, there are orthogonal states in S that are arbitrarily close in their wave functional 
values on continuous A fields. Because of this feature, the space is said to be inseparable. 
Detailed comparisons have been made between Schroedinger and polymer representations in 
the context of quantum mechanics [7] [8] . It is shown that polymer representation, which is 
inseparable in the same way as loop representation, leads to a quantum mechanics that is 
unitarily inequivalent to Schroedinger theory. The cause of the inequivalence is the insepa- 
rability of the Hilbert space obtained from the polymer representation. For the same reason, 
we expect the quantum theory built upon S to be unitarily inequivalent to a theory built 
upon Fock space representation. 

2.2.3 Operators in S and Quantum Geometry 

The natural elementary operators in S are constructed from the flux and holonomy variables, 
respectively defined on oriented surfaces and paths embedded in M. Let S be an oriented 
surface in M, and let f l be an SU{2) valued test function on S. Then a flux operator is 
given by [5] [I] [S] : 

F f {3) = l -f PEte abc dS bc = / fEfdSa = -th8n(G/c 3 h f P^dS a (2.14) 

J S J S J S Ob 

where dS bc denotes the oriented differential surface element; the last equality comes from 
the fact that in the connection representation, E t can be expressed as a functional derivative 
(see (2.3)). Let e be an oriented path in M, and let gl be an SU{2) test function along e. 
Then a holonomy operator is given by [5] [4 J [6J : 

hHe) g ee gliVeitp J _de b A^t^ (2.15) 

where de b denotes the oriented differential path element; in the connection representation a 
holonomy operator is just a multiplicative operator. 
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Specifically, these elementary operators act on a cylindrical function of a single holonomy 
in the following way 

h {j) (e 2 ) g ■ h^ie^A] = h^gh^AA] 

(2.16) 

F^S) ■ hMfaMA] = ^sgn(S,eOf (p)^ 

where sgn(>S, e±) is equal to +1 if e\ intersect S with the same orientation of S at a point p, 
— 1 for the opposite case, and otherwise. The paths e,+ and e 1r result from dividing e\ into 
two new edges from the point p (fig. 1(a)). Recall that spin network states are cylindrical 
functions built from invariant products of holonomies, so (2.15) also determine the actions 
of the operators on S. 

The holonomy operator fo^(e 2 ) fl acts on a spin network state in the straightforward way: 

h^(e 2 ) g ■ S^A} = (^NA^M ( 2 - 17 ) 

The result in (2.18) is a linear combination of spin network states that base on graphs 
covering both 7 and e 2 . Thus the action is in general graph-changing by the inclusion of a 
new edge e 2 . When 7 overlaps with e 2 , the action (2.18) also changes the spin representations 
by coupling M^(e 2 ) g to the existing holonomies in S^j uin . 

As a differential operator, Ff(S) acts on a spin network state following (2.15) and Leibniz 
rules. An important example is given by the case when 7 intersects with S at one and only 
one node v\. In this case we have 

{N v N e ^ 

(g)*n(g)/^) [A] (2.18) 
n i ) 

where i(S, e,/) is equal to +1 or —1 when is above or below S in an infinitesimal neigh- 
borhood of v\. The e n> j/ (with v n e e«/) denotes the oriented path that start from v n and 
overlaps completely with the edge ey (fig.l(6)(c)). The action of Jj(e) with an arbitrary path 
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Figure 1: The upper figure pictures a generic part of an embedded graph, which is an 
embedded edge e\ connecting the two nodes v<i and v 3 . The lower figures demonstrate the 
corresponding definitions of e-, + , e 1r , e 3 1 and e 2 \. In this case, e 2 1 is the same as E\. 

p p ' ' ' 

e on a cylindrical function is given by 
Ue)-^j[A] = ^j>{A] 

f'(h(ei), h(e 2 ), h(e 3 ), h(e Ne ))[A] (2.19) 
= ihGj2[S^ k ■ d*{e k ) - 5^ ■ d^{e k )]f{h{e x ), h(e 2 ), h(e 3 ), h(e Ne ))[A] 

k 



where <9f(efc) and df-(ek) denote the left and right Lie derivatives on h(ek), with the ith 
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generator. Applying it to S^ t j itin one obtains : 



= Inv lihG £ f (g(D)) fee-, • h^\e k )r^ - 5 , ■ r^h^\e k ) 



N„ N K 



*n(g)/l Wl) fe)M4 



(2.20) 

where e(o) and e(i) stand for the origin and the end points of e. 

Combining flux and holonomy operators, one can construct operators that characterize 
the geometry of space. In particular, the spatial area and volume operators have been rig- 
orously constructed in loop quantum gravity. Moreover, their spectra have led to significant 
understanding of the possible Planck-scale structure of space. 

The area operator in loop quantum gravity is derived from regularizing the classical 
expression in terms of the flux variables. The resulting operator for the area of a surface S 

A § \S^ tn ) = y/Fi(S)HS)\S^) (2-21) 

Remarkably, spin network states are eigenstates of A§ for any surface S. Actually, one can 
easily check that 



(2.22) 



where we introduce Planck length through l 2 v = h(G/c 3 ). The summation above goes over 
all points p where 7 intersects transversely with S, and ji represents the spin carried by the 
edge Ci p containing the point p. From above, it is clear that the spectrum of area operators 
in S is discretized, built from multiples of small quanta of order of I 2 . 

The volume operators in loop quantum gravity also result from the regularization of 
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classical expressions using the flux variables, but with much more complicated technicalities. 
Also, there are a variety of different regularization approaches that lead to different operator 
forms for volume [IU]|llj. Here we will follow the one introduced in [11]. The volume 
operators given by this approach act on \Sjyj u i n ) with a set of flux operators over the surfaces 
intersecting 7 only at the nodes. With the help of (2.18), the explicit action boils down to 
EH: 



(2.23) 



Qv n \S7,ji,i n ) — ^ ] Sgn(e n] j, e n j, e n ,k)t Pq ■Jp(&n,i)Jq(&n,j)Jr(&n,k)\S'yJ u i n ) 

where J was defined in (2.19) and (2.20). It has been shown that the spectrum for the 
volume operators are also discretized in S, built from multiples of small quanta of orders of 
lp. By (2.23) the volume operator for a region R acts only on the nodes of a spin network 
state contained in R. Further, each node-wise operation on the spin network state is given 
by the sum of the triplets of the action (2.19) on the edges merging at the node. Therefore, 
the volume operator replaces the intertwiner for each of the nodes in R with a specific linear 
combination of others, preserving both 7 and j{. Consequently, there is an volume and area 
eigenbasis for S, consisting of states with definite 7 and ji. This diagonalizing basis has the 
following interpretation: a certain quantum of volume is assigned to each node of the state, 
and a certain quantum of area is assigned to each edge of the state. Through this basis, S 
provides a concrete description of quantum spatial geometry - the networks carrying Planck 
sized units of areas and volumes, associated with their edges and nodes. 

Finally, notice that the action of an area operator A§ on a spin network state \S>yj u i n ) is 
solely determined by the edges of 7 intersecting transversely with S. Similarly, the action of 
an volume operator Vr on the same state is solely determined by the nodes of 7 contained 
in R. Such information is insensitive to the details about the embeddings of S, R, and 7; 
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only the connectedness and transversality between {ei,v n , S, R} are relevant in the quantum 
geometry given by \S^ ji!in ). 

2.2.4 Knot space K 

The next crucial step is to solve the momentum constraint. Since M g (V) generates spatial 
diffeomorphisms in general relativity [I][2]|3], to solve the quantum momentum constraint is 
to construct a quantum space invariant under spatial diffeomorphisms. 

Set 11 to be a diffeomorphism on the spatial manifold M, which transforms the connection 
fields ll : A — > A'. The action of jx on cylindrical functions is induced naturally by the 
transformation on the connection fields, through the requirement 

^, f [A] = ^ lif/ [A'] 

for all A. From now on we introduce the abbreviation T = (7,ji,i n )- Since T is effectively 
a graph 7 carrying the quantum numbers j% an d i n with its edges and nodes, it is called a 
colored graph. Accordingly, we will denote a spin network state as \Sf). It follows that the 
spatial diffeomorphisms act on S and S* as 

AiSf } = i^r } ; (Sr\v=(S,-.r\ ( 2 - 24 ) 

where liT denotes the colored graph obtained by applying a spatial diffeomorphism 11 to the 
colored graph T. With the inner product (2.11), there is no normalizable state in S that is 
invariant under (2.23) [5j[4j[6j. However, normalizable solutions do exist in the algebraic dual 
space S* [5][1][6], and they can be constructed by the group averaging procedure described 
in the following [5 J [4 J [6J . 

The method of group averaging [22] (23] is a constructive approach to obtain a physical 
Hilbert space for a system that is constrained due to its symmetry. The procedure starts 
from a unconstrained Hilbert space, which is analogous to the off-shell classical phase space. 
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The states in this Hilbert space are in general non-invariant under the transformations of 
the symmetry group. To construct invariant states, the procedure starts with a state in 
the Hilbert space and averages over the transformations of that state by the symmetry 
group. By construction, the results of the averaging will be invariant under the actions 
of the symmetry group. When this procedure can be properly carried out, the resulting 
elements form a physical Hilbert space which naturally inherits the inner product from the 
original Hilbert space [22j[23j. Also, in many cases, the self adjoint operators in the original 
Hilbert space naturally induce physical observables in the physical Hilbert space [22j[23j. 
The method is well understood for the cases of compact Lie groups, while there are many 
unknowns when the group is non-compact, or not a Lie group. Nevertheless, individual 
successful implementations of group averaging do exist for the latter cases. The construction 
of knot space is one remarkable example. 

The symmetry group of concern here is the diffeomorphism group on the spatial manifold 
M, denoted by diffM- Here we will use P to denote 'projectors' that apply group- averaging 
actions to the states in a Hilbert space. Notice that for every 7, there is a subgroup diffM a 
of diffM that leaves the graph invariant, maintaining the set of all the edges and their 
orientations. There is also a subgroup TdiffM,^ of diffM, -y that acts trivially on 7. We define 
the graph symmetry group = diffM, 7/ TdiffM,^, and apply the averaging over Gm,^ to 

any (5p| G S* as 

(Sr\$G M> , = ( 5 rl£ (2.25) 

where N Gm - is the number of elements in Gm,j, which is finite. Subsequently, we further 
apply the averaging over diff m / diffM ,7 , the group that changes the graph, to (Sf\pG M * as 

(Sr\& diffM = (Sf\&G M ^ dijj M / diff Mrl = ^ (Srl^GM^fa (2.26) 

H&diff M / diff M , 7 

The result of {Sf\¥ di jj M is a finite element in S*, since (Sp |P<i^ M |5*f') receive contributions 
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from the finite number of terms from (2.27) where /ry' matches 7 completely. Further, 
(Sf\Pdiff M is diff M invariant by construction. Note that from (2.26) and (2.27), we have 
{S?\Pdiff M = (S 'nr\P : diff M for any /i G diffM- Since the resulting state is characterized by the 
diffeomorphism class [F] of the embedded colored graph F, which is a knot structure, we will 
write (Sf\Pdiff M = ( s [r]l an d C& 11 h a knot state. The inner product (2.12) induces a finite 
inner product for knot states through 

(S[r]hr']> = (s m \Sr,) = (S F |P^ M |S f ,> (2.27) 

With the above inner product, the set of knot states {(srfi|} provides an orthonormal basis 
for knot space K which solves both Gauss and momentum constraints. 

2.2.5 Quantum Geometry in K 

Intuitively, K is just S without detailed embedding information, and each knot state is an 
equivalence class of diffeomorphic spin network states. Therefore, one would expect K to 
inherit the quantum geometry in S that is also ignorant of the embedding details. However, 
the inheritance is not as trivial as it might seem. 

To start with, we study the actions of Ag and Vr on K. Recall that the action of A§ 
on \Sf) takes place at the intersection between S and F. Similarly, the action of Vjj on 
\Sf) happens at the nodes of F contained in R. Based on (2.24), the action of the area and 
volume operators on a knot state (s[f]| is given by their actions on the components {(S^f 1} 
with (j, G diffM- Consider the nodes of /xF contained in R, letting \i runs over all elements 
in diffM- Since the diffeomorphisms move F around freely, the number of nodes of /iF in R 
(in general) changes with /1 . Thus, Vr acts on different number of nodes for the different 
members in {(S^f |}, resulting in a new set of spin network states that are non-diffeomorphic. 
A similar situation also happens to the area operator A§. Briefly speaking, A§ and Vr are 
not operators in K, since their actions on a set of diffeomorphic spin network states would 
result to another set of spin network states that are not diffeomorphic to each other. 
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Clearly, A s and Vr are not operators in K because S and R does not transform under 
diffM and break the symmetry. Therefore, a natural solution would be defining the 'dy- 
namical' surfaces S = S(T) and regions R = R(T), which satisfy S(/jT) = fjf/j,S(T) and 
R(fiT) = fi'fiR(T) for any fi G diff M , with some // G Tdiff M ^. By construction, S'(r) and 
f have the same differential-topological relation as S(fiT) and /if do. R(T ) also enjoys the 
same feature. Then one can define the possible area operator As and volume operator Vr as 

{s [f] \A s = (S T \A m F diffM 

(2.28) 

( s [f]\Vr = (Sr \ V R(r)^diff M 

for any f G [f]. This shows that A s and Vr are now operators in K. Generally, one can 
promote any operator Oq in S* that is sensitive only to the differential-topological relation 
between f2 G M and every f into an operator Oq in K. Accordingly, f2 = Q(f ) satisfy 
Q(fiT) = fi'fiQ(T) for any /i G cJ^m, and some // G Tdiff M ^. 

( Sm \d n = (Sr\d m F diffM (2.29) 

for any f G [f]. From (2.27) and (2.28), S, R and f2 above are equivalently determined by 
the assignments So = ^(fo), Ro = -R(f ) and Cl = f2(f ) given one specific representative 
r from each [f ]. Note that the assignment of S , R and are completely arbitrary 
so far. This arbitrariness indicates the degree of ambiguity that has to be fixed before we 
can meaningfully describe spatial geometry in K . The ambiguity is due to the removal of 
embedding information from knot states. Note that the definitions (2.20) and (2.22) do not 
suffer any ambiguity, since S and R are embedded in a common background M with any 
embedded graph 7, and the intersections between them are thus clearly given. In knot space 
K, the background M is removed and the reference is missing. 

Since S(T), R(T) and 0(f) transform consistently with f, we might consider them as 
genuine dynamical entities. Specifically, any subset of M specified by the values of dynamical 
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scalar fields will have the desired transformation property. This idea of using dynamical fields 
to specify physical regions of spacetime has a long history and explored in abundant works 
[15][16j[18j[17j[19j. We will follow through this idea in this thesis and show that it is indeed 
valid in our context. 

In summary, the quantum geometry of K should be based on the area and volume 
operators of dynamical surfaces S and regions R. Without the background, the relation 
between S, R and every [f ] has to be individually specified. Since the goal of this thesis is 
to examine the semi-classical limit of loop quantum gravity, it is essential to describe locally 
the spatial geometry. To this end, we will need to define S and R in a physically meaningful 
way, so they truly define operators corresponding to variables in general relativity. 

2.2.6 Hamiltonian Constraint and Physical Hilbert Space 

The remaining constraint to be imposed to achieve the physical Hilbert space is the Hamil- 
tonian constraint. Adhering to the quantum geometry of space, the Hamiltonian constraint 
operator H g {N) in loop quantum gravity, with lapse function N, is built with a discrete 
structure. The operator is defined as the limit lim^o Hg(N), while H g {N) is the quantized, 
regularized Hamiltonian constraint acting on S. The full construction procedure is highly 
technical, for details please refer to [J][5J[6J. 

The action of H g (N) on a spin network state \Sf ) is given by a combination of a specific 
set of flux and holonomy operators, whose embedded surfaces and paths intersect with 7 in 
a particular way. The fiducial sizes of the surfaces and paths are e 2 and e, and the expression 
of Hg(N) in terms of the corresponding flux and holonomy operators gives a discretized 
approximation for the classical constraint H g (N), with e indicating the fineness. 

The removal of the regulator is a subtle issue, for lim e _>o Hg(N) does not exist in S. 
The reason is that the embedded surfaces and paths for the flux and holonomy operators in 
Hg(N) shrink with e. Therefore the embedded colored-graphs for the resulting states also 
change with e. Since any two states in S based on distinct embedded colored-graphs have a 
zero inner product, Hg 1 (N)\Sf') will be orthogonal to Hg 2 (N)\Sp) for e 2 < ei, no matter how 



21 



close the two values are. However, the two states HV-(N)\Sf) and Hg 2 (N)\Sf) do belong to 
the same diffeomorphism class, provided e\ and €2 are both sufficiently close to zero. That 
means (spvj | lim^o Hg(N)\Sf) for any Sj does exist, and thus lim^o Hg(N) converges acting 
upon K C S*. 

Unlike the cases of Gauss and momentum constraints, there is a large degree of ambiguity 
in quantizing the Hamiltonian constraint. This is due to the many ways one can regularize 
the same classical constraint using flux and holonomy variables. Moreover, the different 
versions of Hamiltonian constraint operators have been shown to be physically distinct. 
While the quantization ambiguity remains an issue in loop quantum gravity, in this paper 
we will build our model based on the most standard version described in |4||5j|6], and will 
involve only the most general features shared by all the other possibilities. To describe 
the action of lim^Q Hg{N) = H g (N) upon K, we introduce the following notations (fig.2). 
For each embedded graph 7, we label each node with an integer n, and each of the edges 
connected to the node n by an integer pair (n, z)jj For 7, we denote the node n by v^, 
and the oriented path starting from and overlapping exactly with the edge (n, i) by e? ni . 
Since the number of edges connected to a node varies depending on the node, the range of 
% depends on n. Further, to each pair (v^, e^) we assign an oriented path &l ni \ that starts 
from and lies in without covering the other end point. Lastly, to each pair (e^e^-) 
we assign an oriented closed path el . containing the outgoing path el incoming path 
(e7 .j)" 1 , and another smooth bridge path, such that it gives the boundary of an embedded 
surface who has no other intersection with 7- Using these notations and setting k = 8irG/c 3 , 
2 Note that each edge has two labels since it contains two nodes. 
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Figure 2: The left figure pictures an embedded graph 7, with one of its nodes labeled by n = 1 
and the four edges connected to this node labeled by (1, j). The figures in (a) demonstrate 
the corresponding definitions of v^, e^j and e^j, while the figures in (6) demonstrate the 
corresponding definitions of e7 and e 



the action of H g (N) upon is given by (srpi|iy ff (iV)|5'fv) for any srfi and Sfv, with [I] [5] [6] 
(s [f] \H g (N)\Sr,) 
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where h(e) = h^\e) is the holonomy operator with spin 1/2, and the Immirzi parameter 7 
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is set to be any real number |^] Also, the action of H^(N) on any cylindrical function ip^j is 
given by |1] [5] [6] 



= U 



1 

96 
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i,j,k=l 



V 



(2.31) 



where the total volume operator V is defined as 
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(2.32) 



Notice that the embedding details for the assignments of el a and e^ n 4 ^ for any of the 7 
are irrelevant for (2.29), since (srnl does not distinguish deformations within the same diff'M 
class. It is indeed remarkable that lim^o Hg(N) converges upon diffM invariant states. This 
fact allows us to impose the Hamiltonian constraint upon knot space. However, H g (N) 
in general does not preserve K, for the same reason why the operators A§ and V/j fail to 
do so. Recalled that A§ and violate diffM symmetry because they are smeared over S 
and R, which are not diff M invariant. Similarly, (s[f]|ff g (A' r )|S'f/) would not be equal to 
(srf]|^ g (A")|S' M f/) if N(fj,v%) differs from N(v%). On the other hand, the only reason for the 
violation is that N(x) is generally not diffM invariant. In fact, one can check that we have 
(s^lHg^lSf* ) = (s^lHg 11 ^ \ S ^f) , and thus H g {l) is an operator in K. 

The physical Hilbert space for loop quantum gravity should consist of quantum states 
annihilated by the Hamiltonian constraint operator with any lapse function, whose complex- 
ity is demonstrated above. Currently, the search of such states remains a pressing challenge 



3 Notice that part of the quantization ambiguity lies in the possible choices of the spin representation for 
the holonomy operators in (2.28) through (2.30), which we pick to be |. 
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for loop quantum gravity, and is tackled with a great variety of approaches [4J[6J[20J[25J. 
Many proposals are pursued using the Hamiltonian constraint operators with different vari- 
ations of (2.29). Aside from the canonical Hamiltonian approaches , solving the constraint 
using path integral formalism in spin foam models also provides important insights from a 
different point of view [6]. In our model, we will modify the Hamiltonian operator into a 
graph preserving operator in K. The much simplified setting will allow us to apply group 
averaging method, based on certain concrete assumptions, to construct the physical Hilbert 
space for the model. 

2.3 Loop Quantization with Matter Fields 

The previous sections outlined loop quantization of gravitational fields, which is guided 
by the background independence of general relativity. We saw that the procedure leads 
to the knot space K whose elements are drastically different from Fock states in particle 
theories. We encountered difficulty of describing local geometry in K. Tracing the causes of 
these obstacles, we were led to the classical problem of finding local observables in general 
relativity. However, in the quantum theory the problem is further sharpened by the need to 
describe the superposition of background independent states. 

In general relativity, the internal coordinates approach to define local observables has been 
extensively studied pj)] [E] [H] JTTJ [19] . The approach attempts to describe the gravitational 
fields relatively to coordinates given by values of dynamical scalar fields. The scalar fields can 
be a part of the gravitational field or a part of matter fields coupling to the gravitational field. 
In this setting, background independence is respected because both the gravitational field and 
the coordinate fields are dynamical entities in the full theory. Indeed, the approach is mere 
reflection of the daily life reality: we can only measure one physical system by referring to 
another. However, the precise mathematical prescription for the local observables defined by 
internal coordinates is subtle, and often involves the intrinsic properties of spacetime. Among 
the abundant works in this approach, different fields have been used as internal coordinates to 
describe gravity. The most popular method in cosmology introduces four dust-matter scalar 
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fields coupling to the gravitational field, with one of the fields identified as the clock and 
the other three being the spatial coordinates [15J [16J [18J. In this method, the scalar fields' 
degrees of freedom are absorbed to dynamically break diffeomorphism symmetry, resulting 
in gravitational physical degrees of freedom that live on each matter coordinate point given 
by the scalar fields. Consequently, when the matter back reactions are ignored, one expects 
to obtain gauge fixed general relativity from those gravitational local observables. 

In the hope of recovering full general relativity in semi-classical limit of our model, we will 
introduce the dust-matter method in the loop representation. Loop quantization for a system 
of gravity coupling to general matter fields is itself an extensive body of work [26J[27J[13J. 
It has been shown rigorously that loop quantization is applicable to such general systems, 
and the quantization leads to a knot space including the matter degrees of freedom. In the 
loop representation, the matter sector of knot states appears extremely different from Fock 
states in particle theories. Beyond the Planck scale, a much studied speculation is that the 
non-perturbative prescriptions of the knot states would be approximated by perturbative 
quantum field theory on a smooth spacetime background. In our context, we will simply 
focus on the basic setting of knot states including matter fields, and also the matter operators 
acting on such states. 

2.3.1 Classical Theory with Matter Fields 

The system we consider in this thesis consists of the gravitational field, gauge fields, fermion 
fields and scalar fields. The matter sector has the usual interactions mediated by gauge 
fields, and is described by Yang-Mills theory minimally coupled to fermion and scalar fields. 
Moreover, the gravity and matter sectors are fully coupled, so the theory for the whole system 
is diffeomorphism invariant. As mentioned previously, Yang-Mills theory shares many critical 
features with general relativity expressed in Ashtekar formalism. Such features enable a 
single representation of knot space for the whole system. To facilitate loop quantization, the 
matter sector will also employ a triad formalism in the classical theory. 

Aside from diffeomorphism symmetry and local 577(2) symmetry from gravitational cou- 
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pling, the system has additional local gauge symmetry of group Q due to the matter inter- 
actions. In the triad formalism we use half-densitized, Grassmann- valued spinor fields £*(x) 
and their conjugate momenta 7ij(x) to describe the fermions in the system, where i and i are 
respectively the (spin |) SU(2) and Q indices. The real scalar fields and their momenta will 
be described by </>i(x) and -P'(x) . The matter gauge fields will be described by A^(x), where 
i and a are respectively the adjoint Q and spatial coordinate indices. Finally, the magnetic 
fields corresponding to A„(x) will be denoted as B^(x). Written in the terms of the specified 
fields, the Hamiltonian constraint density for the matter sector contains the terms [T5j 



H f (x) = 



E? 



2 v / det£ 



mV/U + D*{t?t£) + ^> T £ + cc. 



P [ P- 1 

HJx) = - p— (x) + - 

1 ' 2«Vdet£ l ' 2 



V det E 



(2.33) 



x 



HymU) = 



g|g r E? + B»B?] (x ) 



Here k = 87rG/c 3 , Q is the Yang-Mills coupling constant, and the covariant differentiation 
D a is taken with respect to the total connection field A a + A a . Also, rjjk are the spatial com- 
ponents of the flat Lorenzian metric with signature (—,+,+, +). The generalized magnetic 
field B* is defined as: 



B!,(x) = ^a bc (d b A^ c - «9 c Ai + f^A^x) 



where f 1 ^ is the structure constant of Q. Including the gravitational part, the full Hamil- 
tonian constraint for the system, with a lapse function N, is given by the sum of all the 
sectors H(N) = H g (N) + H f (N) + H S (N) + H YM {N). Here the matte terms are obtained 
from smearing (2.34) with N. 
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Aside from the Hamiltonian constraint, the SU(2) Gauss constraint and momentum 
constraint also acquire new terms because of the inclusion of the matter fields [27J[26J, 
although we refrain from listing them in this thesis. Further, there is also an additional 
Gauss constraint coming from the gauge symmetry of the group Q. The four constraints 
thus govern the background independent theory of the whole system. 

2.3.2 S and K with Matter Fields 

The classical expressions above suggest that we may generalize the holonomies of SU(2) 
group to those of SU{2) x G to incorporate the gauge fields. On the other hand, since the 
generalized holonomies parallel transport £y(x) and 0i(x), the fermion and scalar fields should 
play roles similar to intertwiners in the loop representation. The loop quantization procedure 
does apply to the classical theory in this way [26J [27J [13J. The result is the generalized knot 
space K that solves all but Hamiltonian constraints. Here we will give a brief description of 
the knot states in K. 

A spin network state for the system is a wave functional of the configuration fields, 

constructed from the loop variables h^\e), h^e), 6^ d \v) and h^(v). Here h^e) denotes 

a Q holonomy in the assigned representation i. The scalar fields are described by h^ k \x) = 

exp(0(x)M k )j), which are called point holonomies in their representation k. The fermion 

fields are described by the irreducible tensors 9^(v), each given by a Grassmann monomial of 

£|(i>) of degree d (for details refer to |T3j ) . Constructed from the above variables, each of the 

spin network states is specified by [26] [27] [13] : 1) an embedded graph[^7; 2) an SU{2) spin 

representation ji and a Q group representation ] i assigned to each edge; 3) two generalized 

Clebsch-Gordan coefficients (intertwiners) i n and i n for SU(2) and Q respectively, another 

Q group representation k n , and a degree d n of the Grassmann monomials assigned to each 

node. Analogous to the pure gravitational case, the wave functional forms of the spin network 

4 Because of the presence of the matter fields, the graphs for the gauge invariant states no longer have to 
be closed or connected. 
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states are defined with cylindrical functions as 



^7,/o i , i „,d n ),a i ,i„,k„) \ A > A, 0, 0] ^m) 

N v N v N v N v N e N e 

Inv { (g) i n (g) i n (g) 0^(*; n ) (g) /i( k ")(t; n ) (g) /^(e,) (g) h^fc) |> [A, A, 9, 0] 

n n n n i i 



where Inv{...} denotes the SU(2) ® ^ invariant contraction. Thus by construction, the 
spin network states are locally SU(2) ® Q invariant. The inner product of these states are 
defined similarly to the previous case, using Haar and Berezin measures for the bosonic and 
fermionic sectors. Using such an inner product, all spin net work states {Sy^j uin: d n ),(u,in,kn)} 
form a (non-orthogonal) basis for the space of the wave functionals satisfying both of the 
Gauss constraints. One can construct an orthonormal basis for S, whose each element is a 
linear combination of the spin network states differ only in fermionic sector. Therefore, an 
orthonormal basis is given by {5 , 7,(j i ,j n ,fe n ),(j 4 ,i n ,k„) = Sr }, obtained from recombining d n into 
k n . Again, a knot state is obtained by group averaging the states in S over diffeomorphism 
transformations 



(Sf\i» d iff = (s^l (2-35) 

The set of knot states {sp]} forms a orthonormal basis for the knot space K for the 
system. By construction, K solves the SU(2) Gauss constraint, Q Gauss constraint, and 
momentum constraint. 



2.3.3 Matter Operators in S and K 

In addition to the gravitational flux and holonomy operators given previously, the opera- 
tors in S also include matter operators that act on the spin network states. Similar to the 
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gravitational operators, the matter operators are also smeared with test functions g, f, g 
and g with appropriate group values. From the gauge fields, we have the holonomy oper- 
ators h \e) g , and the corresponding operators Jf(e) that give the flux of the gauge fields 
analogously to (2.19). From the fermionic monomials of degree one, we have the spinor oper- 
ators 9(p) g = 9^ l \p) g and their conjugate momenta operators fj{p)g = i§^(p)g. Finally, from 
the scalar fields we have the point holonomy operators h^(p) g and the conjugate momenta 
operators pi{p). 

While the physics in the matter sector under loop quantization is an important subject 
discussed in other works (see, for example, [Hj), the goal of this paper is to explore the 
semi-classical limit of the gravity sector, involving the matter sector only to provide coordi- 
nates. Nevertheless, the framework of the matter sector in loop representation is essential for 
obtaining a faithful extraction of the full theory, leaving out only the matter back reactions. 
The important feature of the matter sector for our purpose, is that the fermion and scalar 
operators act on the nodes. 

Four matter scalar fields are required to coordinatize spacetime in the classical theory. 
Moreover, the triads formalism introduces additional SU (2) components of the gravitational 
field that need to be described in terms of matter frames. In the quantum theory, such 
coordinates and frames would be provided by matter operators. Since loop quantum gravity 
is based on flux and holonomy variables, the required coordinates and frame operators are: 
1) four scalar operators for the spacetime coordinates; 2) three SU(2) vector operators and 
their conjugates, for the description of flux variables; 3) two SU(2) spinor operators and 
their conjugates for the description of holonomy variables. 

Clearly, the requirement can be fulfilled provided enough species of fermion and scalar 
fields exist. To obtain the coordinate and frame operators in K, we will define them with a 
generic dynamical point p = p(T) following the notion in section 2.2.5. From the ingredients 
provided by the fundamental matter operators listed above, one can easily construct charge- 
less scalar operators of the form 4>{p), chargeless vector operators of the form J{p) 1 and their 
conjugates, and chargeless spinor operators of the form U(pY and their conjugates. For an 



30 



explicit example, we can set: 

Apt = v}(r k M( P ); J(p) k = f§(n)W(p) ; u(pf = (cRfeiejip); u(p)- k = (cg) s ^t efip) 

Here g 1J is the adjoint-^ Killing-Cartan metric, and (C- 1 'i) k = (CH)^ is the Clebsch-Gordan 
coefficient that projects the direct spinor products into a Q scalar and SU(2) spinor represen- 
tation. It should be clear that there are many other ways to obtain the chargeless scalar and 
current operators, as long as the results have the the desired transformation properties. As 
mentioned in 2.2.5, there is excessive freedom in the definition of p due to diffw symmetry. 
However, as we will describe in detail, the symmetry between different choices of p would be 
dynamically broken by the quantum state of the system through the coordinate and frame 
operators, in an analogy with the classical case. 

Lastly, the matter Hamiltonian constraint operator H m (N) = H f(N) + H S (N) + Hym(N) 
has been also constructed in the loop representation 1 1 3 J . Similar to the construction of 
H g (N), the quantization of H m (N) starts from the regularization of (2.32) using the matter 
loop variables, and results in the well-defined operator H m (N) acting upon K after the 
removal of the regulator. Following the notation used in (2.31)-(2.33) the action of H m (N) 
on a knot state has the form 

(s m \H m (N)\Sr,) = ( S[P] | H^N(vZ)\Sr,) {2 36) 
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where the implicit form of Hm is given by [33] 

( S[f] i^Vf') 

^DP(4m)(1))^^ 

(2.37) 

where e(i) is again the target point of e. We will skip the detailed expression is in this paper. 
Nevertheless, from the operators involved, we know that the action of Hm changes the 
graph and representations (of both SU(2) and Q) of the state. Just as the case of H g (N), 
we have (srfq|iJm"^|Sf ) = (srf'il-ffm^lS^p)- Therefore, H m (N) does not preserve K only 
because N is generally not diffM invariant. 

3 The Graph-Preserving Model 

We start with the Hamiltonian constraint operator for our model, which is modified from 
the standard self- adjoint form 

H(N) = \{H g {N) + H m {N)) + \{H g {N) + H m (N))^ 

where each term is described in (2.30)-(2.32) and (2.36)-(2.37). Recall that the standard 

constraint operator does not preserve K due to the fact that N(x) is not diffM invariant. 

In our model, we will replace the embedded lapse function N(x) with the dynamical lapse 

function N p (pk) defined on a set of dynamical spatial points p = {pk}, where the integer 

k ranges to infinity]^] Using the notion introduced in section 2.2.5, the set of dynamical 

nodes {pk} consists of pk = pk(T) satisfying p k (pT) = f/fxp^T) for any /i G diff M and some 

5 The spatial manifold M contains uncountably many spatial points, but we use only countably infinite 
set {pk} in correspondence to the discrete structure of K . 
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// G Tdiff M ^. Naturally we require that for each f and n, = pk(T) holds for exactly 
one k valuej^] Notice that there are infinitely many distinct sets of dynamical spatial points 
satisfying the above, and we encountered exactly the same ambiguity resulting from the 
background independence in section 2.2.5. Respecting the background independence, we 
will consider all lapse functions N p based on all possible p. The Hamiltonian constraint 
H^^Np) in K is defined by 

{s [f] \H\N p ) = {s [f] 



\{H*{N P ) + H* m (N p )) + \{H;{N P ) + H* m (N p )) 



where we have 

(s [T] \(H* g (N p ) + H* m (N p )) ee (s r | £ + H^)N p ( Pk \ Pk{f)= ^¥ d iff (3.2) 

for any f G [F]. Thus H*(N p ) is now an operator in knot space. It is important to note that 
the classical counterpart H g (N) of H g (N p ) is different from H g (N), which is the counterpart 
of H g (N). In the classical theory, we would expect {H g (N), M g (V)} = 0, since N is now 
dynamical and transforms as a scalar field under the diffeomorphisms generated by M g (V). 
On the other hand, the algebra (2.6) applies to H g (N), with N being a non-dynamical 
Lagrangian multiplier. 

Referring back to (2.30)-(2.32) and (2.36)-(2.37), we note that H*(N P ) is graph changing 
since H g and Hm contains the holonomy operators based on ej ni j\ that is not contained 
in 7. To obtain a graph preserving Hamiltonian constraint operator for the model, we now 
replace e7 a with the full path e^j that coincides with the edge (see fig.2). Also we replace 
^Jnij) w ^h ^Jt,i,ji a dosed oriented path covered by 7, with the outgoing and incoming 
e^j, which contains a minimal number of edges of 7]^] Applying the substitutions to H*(N P ), 

6 Note that this is obviously true had we taken {pfc(f )} for any f to be the set of all points in M. Here 
we define {pfe(r)} to be a countably infinite subset of M, but insist on maintaining the condition. 

7 If such a path is not unique, then i j represent averaging over all the minimal closed paths that satisfy 
the requirements. 
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Figure 3: The lattice rectangular prism 7| C M 3 with the vertices represented by the black 
dots at the locations {X n }, and the links represented by the lines between two adjacent 
vertices. The lattice torus Ttorus for A can be obtained by identifying the opposite boundary 
faces of if. As an illustration for chapter 8, the construction also assigns the location X n of 
a vertex to the corresponding vertex of the lattice torus as its coordinate value. 

one obtains the graph preserving Hamiltonian constraint operator 

H(N P ) = \{H a {N p ) + H m {N p )) + \{H g {N p ) + H m (N p ))* 

which is defined in the subspace Kj- of K based on a single network topology T. 

Consider a lattice rectangular prism in M. 3 with a large number of vertices. Identifying the 
correspondent boundary vertices in the opposite faces of the rectangular prism, one obtains 
a lattice torus Ttorus with N v nodes (fig. 3). By construction, Ttorus has locally cubical lattice 
topology with six neighboring nodes around each of its nodes. With the graph-preserving 
Hamiltonian constraint operator, from now on we will consider only the subspace Kj- torus = A 
for simplicity of the topology. To explicitly write down the action of H(N P ) on A, we will 
re-express the operator using the dynamical nodes and paths instead of the embedded ones. 
Define a set of dynamical nodes {v m = v m (T)}, with m ranges from 1 to N v , satisfying 
v m (T) = v^t-s such that n{m) is one-to-one. Corresponding to the given {v m }, define the 
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dynamical paths {e mJ = e mJ (T)} satisfying e mJ (T) = e^ (m) fc(j) and v m (T) = v^ m) such 
that j(k) is one-to-one. The dynamical closed paths e m ^ = e mt ij(T) may then be defined 
accordingly. Note that in this case the minimal loop e n ^j{Y) covers exactly four edges because 
of the local cubical structure. The action of H g (N p ) on any (spjl £ & can be re-expressed as 



{sn\H g {N p ) 

(Sr\H g{ r)(N p )F d iff 

(Sr\H%fl(N p ) + (^«7)" 5 ^^E^(^lp fc (r)=^(r)) E sgn (e m ,(f ), e mj (f ), e m , fe (f )) 



M*w(r))L, ^5n(i),V( 



i,j,k=l 



x (V(e mj (r)) 



(T) 



Me mj (r)) _, < (F) (1),V( 



(T) 



P, 



(3.3) 



for any T E [r], where we have 



i,j,k=l 



x (/^(e m , M (f )) - h~\e mM {T))) l _ (h-\e m , k {t)))[ ■ \{h{e m , k {Y ))) ' , Vfo 



(3.4) 
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Here, the total volume operator t^f) is defined as 



(T) 



— s S n (e m ,i(r),e mJ (r),e mifc (r)) e P9r J p (e mji (r))J (? (e m)j (r)) J r (e m , fc (r)) 



i,j,k=l 



(3.5) 

Note that the particular choices for {v m } and {e mj -} are irrelevant in (3.3)-(3.5) since they 
serve as dummy variables. The explicit form of H m (N p ) acting on S[f] has a similar structure 
as above, and in brief we have 

(s [f] \H m (N p ) 
= (Sr\H m{F) (N p )P d iff 



Vm 

H { Z m{T)) (^(^(r)), ^(^(r))? , ACe^ij-Cr))? , JiC^cr)), ^(^(r))!, ^(^.^-(r))! 
^ m (r))p(e mjl (r) (1) ^ 

Pi(Um(f))>Pi(em 1 i(f )(!))) Pdi// 



(3.6) 

for any f G [P]. 

Observe that in (3.3)-(3.6), the operators H g( f)(N p ) and H m ^(N p ) act on (Sf | using the 
operators defined with the embedded paths and points that lie inside of f , without disturbing 
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its topology. Consequentially, the states (s^\H g (N p ) and (s^\H m (N p ) will still be based on 
Ttorus and belong to A. Therefore, we can construct our model upon ^ using the self-adjoint 
Hamiltonian constraint given by 

(s m \H(N p ) = (S T \H {f) (N p )F dtff = ( S p,|i(4(JV p ) + H m (N p )) + \{H g {N p ) + H m (N p ))^ 

(3.7) 

The physical Hilbert space for the model will be the space annihilated by H(N p ) for all N p 
based on all p. 

4 Physical Hilbert Space 

The idea of applying group averaging method to solve the Hamiltonian constraint in quantum 
gravity is an ongoing project [24][25j. Specifically, the possibility of solving the standard 
Hamiltonian constraint (2.31)-(2.33) in loop quantum with this method has been broadly 
discussed. While difficulties due to the complexity of diffeomorphism group still await to be 
overcome for the full theory [22] [25], the procedure has been carried out rigorously in the 
minisuperspace models [24][28][29j. The graph-changing operations are frozen out in H(N p ) 
for our model, and such simplicity allows one to solve the constraint with group averaging 
method, based on few concrete assumptions. In the following, I apply this method to the 
graph-preserving model in close analogy to Marolf's approach [28] [29] for the minisuperspace 
models. 

Consider the set of unitary operators {exp(iH (N p ))} given by all lapse functions {N p } 
based on all possible p. To solve the Hamiltonian constraint, we look for an object that is 
invariant under the action of this set of operators. Following the experience in (2.2.4), we 
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look for the group generated by {exp(iH(N p ))}. Explicitly, we consider 

max 1 / 

J] eMiH(N Pk ))\ = {u(g)} g&G (4.1) 

where N Pk for each A; is a member of all lapse functions {N p } based on all possible p, and /= 
means that we identify two expressions if they give the same operator in .ft. G denotes the 
group faithfully represented by the set of unitary operators generated in this manner. Note 
that in our model, exp(iH(N pi )) exp(iH(N P2 )) is not equal to exp(iH(N pi + N P2 )). This 
non-Abelean nature of G is due to the fact that the actions of H(N P ) on two adjacent nodes 
of a knot state do not commute, as a consequence of the graph-preserving modification. 
Assuming the existence of a left and right invariant measure dg for the group G, then the 
desired invariance can be achieved by the group averaging operator over G 

P= / dgU(g) (4.2) 
Jg 

Recall that the group averaging operator Pdiff maps spin network states {\Sf) G 5} into 
knot states {(<Sp|]Pdi// G S*} that span K. In minisuperspace models, there are group 
averaging operators P m ini (corresponding to P in our model) generated by the symmetrically 
reduced Hamiltonian constraint operators. Like P^;//, P m im maps the dual states of the 
states in kinematical Hilbert space into the physical states that lie in the dual space. Here 
we assume that the same applies to our model, that P maps from knot states {(s[f]| £ 
into {P|s[fj) G ^*} that span the physical Hilbert space. 

Since P is invariant under the action of G, so is the element P\tp) G A* with e &■ 
Further, the inner product between any two states = P^i) and |^ 2 ) = PlV^) may be 
defined as 

(*l|*2) = ^l|*2) = ^l|%2) (4.3) 

To be clear, the assumptions being made here are: (1) the existence of the left and right 
invariant measure dg; (2) the operator P maps {(s[f]| G A} into {P|sjf]) e A*}. These are 
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known to hold in minisuperspace models [24] [28J [29J, but remain to be proven in the model. 
Finally, under these assumptions we obtain the physical Hilbert space EI that is spaned by 
{P|srfi}}, which solves the Hamiltonian constraint H(N P ). 

5 Local Observables 

Now we wish to find the local observables in the physical Hilbert space EI that have classical 
interpretations, in order to explore the semi-classical region of the model. In the classical 
case, the values of matter fields can be treated as internal coordinates that enable the 
description of local observables [15] [18] • In the quantum theory, the same idea has been 
applied in different approaches to a varying extent. In particular, a variety of minisuperspace 
quantum cosmology models have successfully reproduced correct low energy limits with 
observables defined by matter clocks (3lJ[30j. 

In this model, we aim to construct the observables that are completely localized in space- 
time in order to make connections to full general relativity. The strategy is to first identified 
spatially local operators in ^ using three scalar matter operators to assign the positions, 
then promote them into localized observables in EI using another matter scalar operator to 
assign the time. Note that the dimension of the semi-classical spacetime is thus put in by 
hand in the model. 

For each set of dynamical nodes {v m } we will construct a set of self-adjoint, commuting 
matter operators, containing the scalar operators {(f> (v m ), 4> l {v m ), 4> 2 (v m ), <fi 3 (v m )}, current 
operators {Vj(v m ), UUv m )} and conjugate current operators {V/ (v m ) , LfJ (v m )} (I = 1,2,3 
for the vector currents; I — 1,2 for the spinor currents). The three scalar operators 
w{v m ), 4> 2 {v m ), 4> 3 (v m )) = 4>{v m ) will be the spatial coordinate operators, while the op- 
erators {Vf(v m ), Uj{v m )} and {V/(v m ), Uj(v m )} will be the spatial frame operators, and 
(j)°{v m ) will serve as the clock. In our model, we use the fermionic sector to construct the 
spatial frame operators, and the bosonic sector to construct the spatial coordinate and clock 
operators. Specifically, we will set the matter gauge group to be a direct product of Ng 
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unitary groups, so that Q = Qn and k n = ( k* , k 2 , k^ 6 ) in (2.35). Then, define 
the spatial coordinate: 

which would be self-adjoint and diagonal in the knot states basis for any real functions 
(0 1 ,0 2 ,0 3 ). Also, assuming there are Nf > 5 species of fermions in the matter sector, we 
denote each species by 9j ( J = 1, 2, Nf). In this setting, we will use the examples given 
in section 2.3.3 for the current operators: 

Jivmt = (vi^mKv^ ; U(v m )j ee idAf^uMu^iv^ 

In classical Hamiltonian general relativity, we obtain diffu invariant local variables by smear- 
ing the field variables with delta functions of the spatial matter coordinates. Setting 0(v m ) 
and 0'{e n> i) to be some gravitational operators in A, we accordingly smear them over all 
the dynamical nodes with a regularized delta function 5 e (<p(v n ) — X) to obtain the quantum 
analogies. Separately, 0(v m ) will be smeared with 5 t ((p(v n ) — X) that picks up the node with 
coordinate value X, and 0'(e nji ) will be smeared with 8 e (4>(v n )—X) and S e ((/)(v n ) — (X+AX)) 
that pick up the path connecting the coordinate values X and X + AX. Here we choose 5 e 
to be a Gaussian distribution with width e, which is finite in correspondence to the discrete 
structure of Explicitly, for any S[f] G M. 

(s [T] \6(X) = (5 r |^6(t; n (f))det(A0(t;„(f))) ( 5 e (0(t;„(f)) - X)F Mff 

n 

(smlO'(ew) = (^rl ^0(e n ,(f )) det(A0(^ n (f ))) det(A0(t; (n , ) (f ))) 

n,i (5.1) 

x^(4>(v n (T )) - X)^(0(^)(f )) - (X + AX))P di// 
where the dynamical nodes {i>( n ,i)} are defined by ■U( ri)i )(r) ee e nji (r)( 1 ) with e(i) being the 
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target point of e. In correspondence to the classical expressions, where the smearing is 
performed through integrations over a coordinate space, the discretized version (5.1) of 
smearing is summed over the small coordinate volume elements. The coordinate volume 
element operators are given by 

(^,/|det(A0(^(f))) 

_ _ A A A (5.2) 

= m,f\ 2^ sgn(e nii (r),e ri)i (r),e n , fc (r)) det(A0 e - ni( r), A0 e - nij(r) ,A0 e - n fc(r) ) 

Here, the coordinate difference operator A(f> Sni (T) is given by 

(^7,/l A 0e- nil (r) = (^7,/l[^(n,i)(f)) -^(«n(f))] 

It is crucial to note that the spatially localized operators obtained in (5.1) do not depend on 
the choices of {v m } and {e ni } which serve as dummy variables. Instead, the operators are 
localized by the spatial matter coordinates. 

The classical gravitational fields in Ashtekar formalism are SU(2) tensors, and their 
scalar products with appropriate SU (2) matter currents provide the invariant components. 
In other words, the physical components of the gravitational fields can be described with 
respect to the physical frames of the matter currents. In the model, we use {V}(v m ), U\(v m )} 
and {V/(v m ), Uj(v m )} to extract the physical components. Explicitly, for any S[f] G A 

(5.3) 

(S[f#(en lfe )5- = (5rl^(^(n )fe )(r))^(e n>fc (f))^(^(f))^ // 

Using (5.1), we obtain the spatially localized operators as J(ex,Ax)i and h(ex,Ax)j- Notice 
that while these operators are well defined in A, how well they describe the physics depends 
on how well the matter coordinates behave. 

The next step is to promote the spatially localized operators into fully localized observ- 
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ables in EI by using the clock field operator. Each state in EI represents a quantum spacetime, 
and each state in ^ gives a spatial quantum geometry. Therefore, the temporal localization in 
our model is equivalent to a foliation of a quantum spacetime |f) GH into the spatial slices 
{(■0t| £ labeled by the clock time T. Correspondingly, a localized observable 0(X,T) 
will first project a state G EI into a state \i[)t) = IItI^) £ which is the dual state of 
(■0r| G Then, it will act upon j^r) by the spatially localized operator 0(X). Finally, it 
will apply the group averaging operator P to bring the result back into EI. With the clock 
operator 4>°{y n ) that has continuous spectrum in A, the localized observables described above 
have the explicit forms: 

6(X,T) = PO(X)fl T (5.4) 

where fir is the operator that projects onto the eigenstates of (j)°(v n ) with the eigenvalues 
lying within a range e around T. In our model, the action of IIt on ^ G EI is given by 

( S[ r]|nr|*) = (Sr|n T( r)P*//|*> = (5r|^m|n^o(^(f))5 e (0 o (^(f)) - T) J P di// |tt) 

(5.5) 

for any s^p] e & an d f G [f], where u^o(v n ) is the Hamiltonian speed of the clock 

(5.6) 

Also, since in general the operators z>^o(v n (r)) with different n do not commute, we fix the 
ordering ambiguity by applying the symmetrization sym{...} in the ordering of n. Clearly, 
the choice of {v n } is again irrelevant. The classical counterpart of P^oiyniY)) is i^o(x) = 
{4>°(x),H(N = 1)}. When all fields are on-shell, H(N = 1) generates diffeomorphisms 
given by unit normal flows orthogonal to spatial slices, which transform 0°(x) = 0°(x, 0) 
into (f)°(x,t). Therefore, when all the constraints are satisfied, z^o(x) is simply the speed 
d t (p (x, t)\ t= o of the clock with respect to the proper time carried by the flow. 

Finally, the localized gravitational observables are constructed by promoting J (ex, ax) i 
and h(e x ,Ax)j using (5.4) into J(e x ,Ax,T)j and h(e x ,Ax,TYj. 



(SrM^f)) = (S r \- i°(v n (r)),H {F) ( 
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6 Matter Coordinate Conditions 

In order to give sensible descriptions of the physics, the localized observables must refer to 
reasonably well-behaved internal coordinates. This requirement is to be met by the choices 
of the scalar fields, current frames, and the quantum state together. 

6.1 Conditions on the Clock 

As mentioned in the previous chapter, the purpose of the clock field 0° is to physically define 
a slicing for a quantum spacetime. To explore the dynamics in the model, 0° has to provide 
a slicing that reveals the causal evolution of the system. In classical general relativity, such 
a clock specifies an ADM foliation, in which each slice provides the full information of the 
whole spacetime. In other words, the field values on a slice at any clock time can be used to 
reconstruct the whole spacetime. 

In the setting of our model, the two spatial slices given by two distinct clock times 7\ and 
T 2 are described by f[ Tl |\l/) and ffr 2 |\[ / ). Analogous to the condition above in the classical 
case, the two spatial slices defined by the ideal clock 0° should contain the same information 
about the space time up to some quantum fluctuations. In other words, the group 
average operator should be able to recover the spacetime from either of the slices 

M Tl |# ) + 0(h) = pn T2 |^) + 0(h) = |^) (6.1) 

Notice that the condition involves not only 0° , but also This agrees with the classical 
theory in which the proper clocks are chosen according to the state of the whole system. 
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Using the definition (5.4), the condition (6.1) implies 

Fd 1 (X)d 2 (X)fl Tl \^) 

= pdi(x)d 2 (x)n Tl pn Tl |^> + o(h) 

(6.2) 

= p6i(x)n Tl pd 2 (x)n Tl |^) + o(h) 

= 6 1 (X,T 1 )d 2 (X,T 1 )\*)+0(h) 

More generally, denoting $(Oi, 2 , On) as an ordered function of N linear operators, we 
have 

p$(Oipo, 6 2 (x), 6 N (x))fi Tl I*)] 

(6.3) 

= Ti), 6 2 (X, 7i), CV(X, 7i))|*) + 0(h) 

In the context of this thesis, we will need to consider the interactions applied to the clock. 
For later use in describing this interaction, we introduce the operator $i n t(v n ): 

M&intM = (Sr\^(v n (r)) [^ (v n (r)),H {T) (l)\ F diff (6.4) 

where u,po(v n (T)) and ^^(^(f )) are the Hamiltonian speed defined in (5.6) and its inverse. 
Note that $j„ ( (t)„) contains mainly the 'acceleration' of the clock, which is given by the 
commutator in (6.4). 



6.2 Conditions on the Rulers and Frames 

In addition to a good clock, we need good spatial coordinates and frames to physically spec- 
ify positions and directions at a given clock time. In classical general relativity, ideal spatial 
coordinates for a certain foliation should provide local diffeomorphisms between M 3 and the 
space. In other words, the coordinatization should be locally one to one and preserve smooth- 
ness and adjacentness. In loop quantum gravity, the states of spatial slices are described by 
knot states with discrete structure. Therefore, in the model the matter spatial coordinates 
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should locally provide a 'discretized diffeomorphism' between the quantum space and a set 
of coordinate values. 

Before imposing the spatial coordinate conditions in the model, we need to show how 
the coordinatization can occur in the first place. To start with, we observe that by applying 
to IItiI^) all arbitrary functions &(J(ex,Ax)i, h(ex,Ax)j, w(ex,Ax)j) of the set of spatially 
localized operators, one obtains a set of states constituting a linear subspace of £*. This 
subspace {$( J(e x ,Ax)i, h(e x ,Ax) j, W(e x ,Ax) |^)} = wm ^e called the relevant 

subspace for fl Tl \fy). Recall that there are a great many distinct sets of dynamical nodes 
and dynamical edges compatible with them. This ambiguity is a manifestation of difju 
symmetry in the quantum theory, which involves the superpositions of the knot states that 
are defined background independently. Here we impose a requirement on the state |^) such 
that the corresponding 80% T dynamically break the symmetry, picking up a preferred {v*} 
and {e* j} by being an approximate eigenspace of {4>(v* L )} and {4>(V( ni ))}, where we again 
define v* ni ^(r) = e* i (r)( 1 ). In other words, we require to satisfy 

det(A0«))^(0«) - X m )\i>) = 5 m ,M + 0(h) 

det(A0«)) det(A0(^ ni4) ))^(0«) - X)5^(v* M ) - X m - AX mj )|^> = <W^#> + 0(h) 

(6.5) 

for any \ip) e My Tl - The appearance of det(A0(u*)) is needed for the same reason explained 
after (5.1). When this is satisfied, the X m is assigned to v* m and the coordinate gap AX m j 
is assigned to e* m ^. Since 0(t>*) is diagonal in the orthonormal basis of knot states, the 
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condition (6.5) is equivalent to 

( S[f] |^> det(A0(O)m<) -X m )\s [f] ) = (s m \1>)6 m „\s m )+0{h) 

(s [t] \i>) det(A0(^))^(0«) - XMkvUi)) ~ (X m + AX m ,))\s [f] ) = ( S[f] |^ m ,„<yV>) 

+0(h) 

(6.6) 

for any dual knot state \s^} and any \ip) e According to (6.6), the condition (6.5) 

is also equivalent to the requirement that 

(s[r]|V>) k (v n )\s [f] ) = (s [f] \i)) k^^lspj) ; 

(V( k^«)),0 2 ( k^«)),0 3 ( kft«))) =X n + 0(h) (6.7) 

(>( k^ n>i) )),0 2 ( k^(^ 4) )),0 3 ( kft^))) =X n + AX n4 + 0(h) 

for any dual knot state |srri) and any e . From (6.7) we immediately see that 
{ v ( m ,i)} = i v n} and {(X m + AI mi ,)} = {X„}. The physical meaning of (6.5)-(6.7) is cru- 
cial: since the matter spatial coordinates are what we refer to in spatial measurements, the 
preferred dynamical nodes {v* } represent the physical spatial points we observe. Consequen- 
tially, there is also a set of preferred dynamical spatial points p* = {p^} agreeing with {u*}, 
such that p* m (T ) = v^(T) for m < N v . It should be emphasized that {v *} and p* are induced 
by ^Ti> so they depend on the spatial quantum state fl Tl |\l/). For this reason, as one can 
see in (6.6) that t>*(r) can be arbitrary for |srfi) that is orthogonal to . 

The spatial coordinate condition can now be imposed on the map (v * , e* J — > (X n , AX Uti ). 

8 Condition (6.6) implies that the coordinate gaps between adjacent nodes are comparable to the width e 
of 5 e , such that there is only one dynamical node corresponding to each coordinate value. This set-up is only 
for simplicity. Coarse grained versions of (6.6) can be imposed when e is much bigger than the coordinate 
gaps, without effecting the following arguments about the semi-classical behavior of the model. 
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Recall that .ft is based on the lattice torus Ttorus- Analogous to the coordinate maps for a 
torus manifold, the map (i>*,e* J — > (X n , AX nti ) in large scale should appear to be smooth 
around most of the physical spatial points. On the other hand, since there is no global 
diffeomorphism between M 3 and a 3-torus, we expect the map to appear discontinuous at 
some specific set of the points {v*} C {v*}. To approximate local diffeomorphisms between 
the physical spatial points and IR 3 , the model imposes the spatial coordinate conditions 
as the following. For every e* A with e* ni (Y) fl {v* b (T)} = 0, we demand (AX^I < d, 
where d represents a small coordinate gap, such that the map appears continuous in large 
scales everywhere except {v*}- For any m the set {AX m j}\ s * jf)nfv* (T)\=d should define 
a parallelepiped in M. 3 up to an error of 0(d), such that the map appears smooth in large 
scale everywhere except {v* nb }. Notice that once (6.7) is satisfied the coordinate conditions 
can be achieved easily through a suitable redefinition of the scalar functions (0 1 ,0 2 ,0 3 ). 

In (5.3), we defined the physical SU(2) components of the tensorial observables with 
respect to the matter current frames. Obviously, valid frames for this purpose should be 
physically orthonormal to each other. Explicitly, we require that for any \ip) G 

(djvivDixM) = 6\{XM) + o(h) 

(d'luju-;)(x n m = d'l(xM) + o(h) (6.8) 

(e^V/V/V k K )(XM) = e IJK \i>) + 0(h) 

for any gravitational operators Oi(v n ) and 0'l(v n ). Matter frames that satisfy these condi- 
tions are orthonormal in terms of the spatial geometry given by the gravitational sector of 
the state \ip). 

With the help of the appropriate matter spatial coordinates and frames, any operators 
acting on the space ^ Tl can be re-expressed in terms of the spatially local operators. It is 
through these new expressions we will decipher the physical meanings of the abstract knot 
states. From now on, we denote X and (X, AX) to be variables taking values from {X n } 
and {(X n , AX nA )} respectively. 
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To start with, we assign physical lapse functions referring to the spatial coordinates. For 
every function Af(X), there is a correspondent lapse function N$ with N^(p* m ) = J\f(X m ) 
for m < N v . Therefore, given any function Af(X) , one can mimic the Hamiltonian constraint 
operator with the lapse function using the spatially local operators. 



H g (M) 



U*{N) - (ihKry)- 
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X] A/ '( X ) X sgn ( e ^- A ^' e ^,AV , ex,A^) 



AX,Ar,AZ 



x (fc- 1 (e w )) r [(Mex,Ax))_ , [nf(l),v] (h'\e XAY )) _ (h(e x , AY )) _ , [wf(l),\>' 
^(ex.Az))^ [(Mex,Az))' V 



(6.9) 



where %f (AT) is 



«^«7) ^X^W X s g n ( e ^Ax,ex,A^> e x,Az) (Me^,Ax,Ay) - ft ^x.ax.a^)) - 
x (/r 1 (e x , AZ ))j'- (^(e X) Az))^V 



and 



1 
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(6.10) 



^ sgn (e x ,Ax, e x ,AY, ex A z) e PQR J (e XAX ) P J {e XAY ) Q J(e XAZ )R 

AX,AY,AZ 



(6.11) 
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Clearly, the operators H g (N^{) and H g (Af) are different operators in the full space .ft. How- 
ever, thank to the coordinate and frame conditions (6.5) and (6.8), for any \ip) e ^ ^ we 
have 



H g (N$)W=H g (Af)\tl>)+0(h) 



(6.12) 



Therefore, H g {N p %) and % g {M) act equally on any IT^ |\&) e ^ ^ up to quantum fluctuation. 

Set a(e YAY ) to be the operator obtained from localizing a(e n) j), which is defined by the 
following. Each knot states (srfi| € .ft is an eigenstate of a(e nji ) with an eigenvalue +1 or 
— 1. The knot state is assigned the eigenvalue +1 if the embedded edge of f overlapping 
with e nj j(r) has the same orientation as e n ^(T), and the eigenvalue —1 if the orientations 
are opposite. When acting on any \ip) e ^, Tl , the spatially local operators also satisfy the 
equations: 



AX I 



J(e X: 



AX I 



+ 0{h) 



J(e x ,Ax)i, h(e Y ,AY)j 



5x,y5Ax,Av^p7( r /)j^(ey,Ay)ilV') 
-Sx,Y+AYS-^ XA Yill'r(T I ) I R h(eY,AY)j\'^) + 0{l 2 p h) 



J(ex,Ax)i, h f (e Y ,AY)j 



H e X,Ax)l,k e YAYYj 



= Sx,Y^AXAY^pl{ri)jh\e Y ,AYYK\^) 

-5x,y+Ar^-Ax,Ay^7(r;)^/i t (ey,Ay)j |V>> + 0{l 2 p h) 

= + o(i 2 p h) 



(6.13) 



h(e x ,Ax)z, h\e Y ,AY)j 



o + o{ilh) 



J(ex,Ax)ii J{ c y,ay)j 



dx,Y S Ax,AYilpl£iJ J(e YAY ) K a(e YAY )\^) + 0{l p h) 
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Again, it should be emphasized that the equations in (6.13) are not equations of the 
operators alone. The equalities are derived from the conditions (6.6) and (6.7) satisfied by 
the state flrj^), so they arise from the physical conditions of the system. 

7 Coherent States 

We have required the state to satisfy quantum coordinate conditions, such that the set 
of local observables {J(ex,AX)^);} and {h(e x ,Ax, T) j} are expected to give a meaningful 
description on the gravitational sector, with T around the moment Ti . Since our goal is to 
obtain the semi-classical limit, specific conditions of coherence should be imposed on |\&). 
As in the case of quantum cosmology, we need semi-classical conditions on the gravitational 
local observables and the momentum of the clock. Therefore, the state |\&) is required to 
satisfy 

J(ex,Ax,Ti)i\V) = (<*|J(ex,Ax,7i) 7 |*)) |tf> + 0(h) 
h(e x ,Ax,Ti)j\iff) = ({*\h(ex,Ax,T 1 ) I j\*)) \*) + 0(h) (7.1) 
^(X,7\)|*> = ((^(x,^)^)) \*)+o(h) 

Additionally, we also expect the expectation values to appear continuous in terms of the 
spatial coordinates for a semi-classical state. Therefore, for any two e* j(F) and -(r) that 
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share a common node and form a smooth path, we will impose: 

miexM^Yji*) = mkex^Mim+oid) (7.2) 

(*|^(X n ,71)|*) = (*|i> p(X m ,T 1 )|*) + O(d) 

In this thesis, we will assume the existence of the states that satisfy (7.1) and (7.2) without 
explicitly constructing them. Because of the appearance of K or l 2 v in every result of (6.13), 
this assumption is rather mild. Note that the coherence conditions are defined with respect 
to the observer who carries the clock 0°. The conditions say that the quantum spacetime 
has a sharply defined spatial geometry and extrinsic curvature on the spatial slice at 
clock time T\. Therefore, the quantum spacetime is expected to be semi-classical around 
that moment. 

8 Emergent Gravitational Fields 

To make contacts with classical general relativity, the expectation values in (7.1) need to 
be interpreted in terms of classical fields. The matter coordinates provide a natural way 
to make such a translation from the discrete structure of knot states into the continuous 
configuration of classical fields. 

For an explicit example, we will first pick a simple spatial matter coordinate system. 
Recall that the lattice torus %orus can be constructed by identifying the opposite boundary 
faces of a lattice rectangular prism J| C IR 3 , which consists of the vertices {V n = X n } and 
links {k}. Here the bars indicate the embedding in IR 3 . Such construction naturally gives a 
coordinate map i>* — >■ V n — X n that satisfies the spatial coordinate conditions (fig. 3). Also, 
the approximated spatial coordinate space is naturally a rectangle region P D 7| inside of 
IR 3 . In this case, the set {v*} where the map is discontinuous is the set that is mapped into 
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the boundary of P. 

Once a matter spatial coordinates are chosen, our model identifies every e* i satisfying 
e* j(f) n {^.(r)} = with an embedded path in M 3 under the guidance of the matter 
coordinate values. Specifically under our choice of coordinates, each such e* A is identified 
with the oriented path ex„,Ax n:i that goes from the vertex X n to the vertex X n + AX nji , 
which overlaps completely with the link connecting the two vertices. Subsequently, the 
model choose a cubical cell decomposition dual to 7|, dividing J 3 into a set of cubical cells 
{cx n }- Each cell cx„ uniquely contains a vertex X n , and the boundaries of all the cells 
consist of a set of smooth faces {si} whose each element Si> intersects transversely with an 
unique link among {1$}. Suppose li> links X n and X n + AX n ^, and denote Sx n ,Ax ni C P 
to be the oriented surface overlapping completely with Sj/ that is dual to and has the 
same orientation as ex n ,Ax ni - The picture motivates an interpretation of the expectation 
values in (7.1) as being given by smearing a smooth gravitational fields defined in J 3 over the 
corresponding elements from {ex, ax} an d {Sx,ax}- Explicitly, we pick a fitting algorithm 
that maps the expectation values {(J(ex,Ax, T)i), (h(ex,AX,T)j)} and {(v^X, T))} to the 
values of the smooth fields {Ef (X, T) , A J b (X, T) } and {z/^(X,T)} defined in P. The model 
requires the fitting algorithm to obey the following rulesj^] 

^o(X,T) = (^o(X,T)) 



/ E^dsa^iJiexAx^h) (8 - L) 

J Sx A 



Pexp[ / Ai{T){rj)de h \l = (h(e XAX ,T))\ 



Due to the lattice torus topology, the continuous condition (7.2) would imply boundary 
conditions on the emergent fields in the coordinate space P. Denote the three pairs of 



9 Note that such an algorithm is guaranteed to exist, since we are fitting the smooth fields with infinite 
degrees of freedom to the finitely many data points given by the expectation values of the local observables. 
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boundary faces of I 3 as (dl~£,dl x ), (dly,dl y ), and (dl^,dl z ), the continuous condition 
(7.2) on the lattice torus implies: 

u 4P (X,T 1 )\ 0I + tVtS = iy (/) o(X,T 1 )\ dI -^ z + 0(d) 

EftX^)^ = E^X, Tl )\ dI - yz + 0(d) ^ 

■ = A J b (X,T 1 )\ dI -^+0(d) 

The choices of J| C M 3 , the cell decomposition {<Sj}, and the fitting algorithm described 
above are restricted but non-unique. However, any chosen set gives a valid correspondence 
between and the smooth fields. 

It is important to check that the emergent gravitational fields transform correctly un- 
der changes of spatial coordinates and frames. First, we consider two observers using the 
same clock and current frames, but two different spatial coordinates given by 4>(v^) and 
0'(i>*)P^| Suppose that at T 1; the two spatial coordinate systems v* — > X n e 7| and 
v * — > X' n G I' z satisfy the spatial coordinate conditions, and are approximately continu- 
ous except at {v*} and {v*,}. For the coordinate gaps, we arbitrarily set d > d'. Any 
e* ni satisfying e* j(f ) PI {u£ (f)} = and e* A (T) fl {«* /(F)} = would be mapped into an 
embedded path in each of the two coordinate spaces I 3 and I' 3 . The two embedded paths 
in the two coordinate spaces are respectively ex n ,Ax ni and ex> n ,AX' .• Also, the cell decom- 
position of I 3 dual to If, and of I' 3 dual I'\ induce surfaces {Sx n ,Ax ni } and {Sx^,ax' } 
in the two coordinate spaces. By construction, the transformation (ex n ,Ax ni , Sx n ,Ax n:i ) ~* 
(ey/ ax' i Sx> ax' ) is contravariant under the coordinate transformation X — > X' up to an 
error of order d. On the other hand, we have (J(ex n ,AX ni i = (J(e-x' n ,AX' ■■> + 0(h) 

10 It would be interesting to consider the case where the two spatial coordinate operators are based on 
different preferred dynamical nodes. In such a case the two observers has different notions of physical spatial 
points. Here we would focus on the simpler case with their agreement on {i>*}. 
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and (h(e Xn ,Ax n , i ,T 1 )f) = (h(e x ^,AX' n 4 , T i)i ) + Oifi). Thus from (8.1), we conclude that the 
emergent gravitational fields transform as they do in the classical theory under a change of 
matter spatial coordinates, up to an error which is suppressed when both h and d are small. 

Next, we consider two observers using the same spatial coordinates but two different 
frames, each satisfying the conditions (6.8). At the clock time T\, the transformation matrices 
between the two frames are given by 

(¥ | V! V'/,(X,Ti)\^) =TZ I j,(X : T 1 ) + 0(h) 

(*\v' i I 'v}(x,T 1 )m = ny\x,T x ) + o(K) 

(8-3) 

mu-lu'jMx^m ^u^x^ + oih) 
Wu'/'ulix^)^) = Uf (X,71) + o(h) 

It follows from (6.8) and (8.3) that K' 1 = 70 and U' 1 = W. Therefore the clock time T\, the 
observables' expectation values for the two observers are related by 577(2) transformations 

^,(X,T 1 )(J( 1 S x ,ax,T 1 ) / ) = (JiSx^T^r) +0(h) + 0(d 3 ) 

(8.4) 

uf (X + AX,T 1 )(/i(e x , AX ,T 1 )f)W|,(X,T 1 ) = (h(e x>AX , T x )f, + 0(h) + 0(d 2 ) 

Again, combining (8.4) with (8.1), we conclude that the emergent gravitational fields trans- 
form as they do in the classical theory under a change of matter frames, up to an error which 
is suppressed when both K and d are small. 

Altogether, we have the transformation between the two emergent gravitational fields, 
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given by two different sets of matter coordinates and frames, of the form 




dX 




(X,T 1 ) + 0(h) + 0(d) 



Vdet E> 



dX a 



VdetE 



(8.5) 

r)X a - 
{AAX'^f; = (tfAMiX,^) +itfd a U(X,T 1 ))% + 0(h) + 0(d) 



which has the correct semi-classical limit. 

The emergent gravitational fields describe the large scale limit of the physical quantum 
geometry we have been searching for in section 2.2.5. Corresponding to each set of {e nji } 
and {v n }, we construct a set of dynamical surfaces {S n> i} and regions {R n } satisfying the 
following conditions. For any f , the embedded surface S n>i (r) transversely intersects once 
with only e n ^{T) among {e m j(T)}, in the same orientation. The embedded region R n (T) 
contains only v n (T) among {v m (T)}. Using (5.1) and (5.4), one obtains the local area 
and volume operators {A(Sx,ax, T)} and {V(Rx, T)}. Moreover, the expectation values 
of the localized geometric observables {A(Sx,ax, T±)} and {V(R x ,Ti)} measures the areas 
and volumes of {S*^} and {-R*} corresponding to {e* ni } and {f*}. The sets {S*^} and 
{Rn} represent the sets of physical surfaces and regions for the state |^) at 7\ under the 
matter spatial coordinates. Further, the coordinate map v* — > V n — X n naturally identify 
the physical surfaces and regions with {Sx n ,Ax nii } and {cx„} in the coordinate space P. 
Finally, the quantum geometry given by the action of {A(Sx,ax,Ti)} and {V(Rx, Ti)} on 
|\&) is coherent, and has a large scale limit conforming with the emergent gravitational fields 




XM^+Oih) 



(8.6) 
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where the left hand sides are exactly the classical expressions of the area and volume in 
terms of the emergent densitized triad fields. 

9 Emergent Constraints and Diffeomorphism Algebra 

Now we move on to explore the consequences of the symmetry of |\&), or equivalently the 
symmetry of P. In particular our goal is to see, in large scales and up to the matter back 
reactions, what equations are imposed by the symmetry on the emergent gravitational fields. 
The right invariance of P implies that the physical state satisfies 

(^|Pexp (iif(eWp)) Aril*) = (tf|M Tl |tf) (9.1) 

where N p and the real number e are arbitrary. Denoting O' =0 + 0^ and taking e to be 
small, we have 

(*|M(Agn Tl |*> = (*|M;(^ p )n Tl |*) + (*|M^(^ p )n Tl |*) = o (9.2) 

Specifically, we can set the lapse function to be N p = N^l corresponding to an arbitrary 
Af(X). Then it follows from (6.12) that 

(*|P7}-(A^)n Tl |*> + (*|M^(7V^)n Ti |^) 
= (*|M-(^)n Tl |*) + (*\FH- m (N«)TL Tl \*) + 0(h) 

(9.3) 

= {^\FH- g (Af)tl T M+0(h) + e m 
= 

for any Af(X), where we denote the contribution from the matter sector as e m . Recall 
that Hg(J\f) is an ordered series of spatially local operators. To indicate this we write 
i-L g {J\f) = H g (J\f, J, h), where J and h denote the collection of all the spatially local operators. 
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Moreover, the condition (6.3) for the clock gives 

= (*\PH- g {Af, j,h)n Tl \V) + e ri 



(9.4) 

o 



The coherence conditions (7.1) says that the local observables in (9.4) can be replaced by 
their expectation values, with errors of orders of h: 

(p\H- a (M, J(Ti), h(Ti))m + 0(h) + e m = mU g {U, (J(rO), (h(Ti)))m + 0(h) + e m = 

(9.5) 

Under the correspondence (8.1), we now rewrite (9.5) in terms of the emergent gravitational 
fields. 

(m- 9 W, (J(T\)>, (h(Ti)})\^) + 0(h) + e m = H g (^)\ Ef{Ti)Ai(Ti) + 0(d 4 ) + 0(h) + e m = 

(9.6) 

where H g (H) is exactly the classical Hamiltonian constraint given by (2.4), written in terms 
of the emergent gravitational fields and matter coordinates 



It follows that up to the corrections coming from the discretization of space, quantum fluctua- 
tions, and matter back reactions, the emergent gravitational fields of |\&) satisfy gravitational 
Hamiltonian constraint. 
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We now move on to show that Gauss and momentum constraints also emerge correctly 
in the semi-classical limit. The invariance of P again implies 

(*|Pexp (iHieNpf) exp (#(ie'M p )) exp (if(-ieN p )) n Tl |*) = (^|Pn Tl |*) (9.7) 

Setting e and e' to be small, we have 

(*\F[H(M p ),H(N p )}n Tl \V) = (9.8) 

for any M p and N p . Setting the lapse functions to be and Nfi with arbitrary M.(X) 
and Af(X), we have 

^\F l -[H(M^),H(N^)]U T M 

= (*|P^ft(M^), J ff-(^)]n 31 |*> + 4 

= <*|p^[w-(jW),w-W]nT 1 |*> + o(ft) + 4 (9-9) 

= 



Here we have singled out the gravity-gravity term of the commutator and denoted the terms 
involving H m as e' m . The commutator \H' g {M), %' g {J\f)\ is then carried out using (6.13): 

(*|p^[?i-(M),?i-(^)]ri Tl |*) + 0(h) + e' m 

= ($?\PC g (M,Af, J, h)IL Tl \V) + (V\W a , a {M,N, 1 k a)fl Tl + O(ft) + e' m 

= mc 9 (M,N, At,), fe(7i))|*) + (*|c^(jw,^, i(Ti), MTi), am))!*) + O(fi) + 4 

= 

(9.10) 



where we separate the result into two terms according to the presence of a. Among the two 
terms involving C g (M) = C g (M,M,J,h) and C g<a (M,J\f) = C g>a (M,Af,J,h,a), only the 
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latter contains a. Similar to the previous case, the semi-classical limit of (9.10) in terms of 
the emergent gravitational fields is obtained as 

(*\C g (M,M, (.7(71)), (MTO))!*) + W 9 , a (M,Af, (AT,)), (MTO), (a(Ti))m + 0(h) + 4 
= {H g (M),H g (Af)}\ EnTi)Ai{Ti) + 0(rf 4 ) + 0(h) + e' m 

= 

(9.11) 

Note that the expectation of a for \^/) is of order of one, since its spectrum contains only 
+1 and —1. Given that, the term C g>a only contributes as 0(d i ). Thus the leading order 
contribution only involves the expectation values that can be translated into the emergent 
gravitational fields. Since M. and Af are arbitrary, it follows from (2.6) that the emer- 
gent gravitational fields also satisfy gravitational Gauss and momentum constraints, up to 
the corrections due to the discretization of space, quantum fluctuations, and matter back 
reaction. 

The above calculations suggest that the operator V, - g (Af) is the quantum counterpart of 
H g (M). According to the classical case, W, g is then expected to generate the translations on 
only the gravitational fields, in the direction perpendicular to the spatial slices specified by 
the clock. Moreover, the calculations also suggest that C g (Ai, Af) is the quantum counterpart 
of {H g (M), H g (M)}, which generates combinations of spatial diffeomorphisms and local 
SU (2) transformations on only the gravitational fields, with respect to the matter spatial 
coordinates and frames. As a first step to validate this interpretation, we list the results of 
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the semi-classical limits of the commutators 



{9\F-[Hg(M),H- g (Af)\iL Tl \v) 



{H g (M),H g (Af)} \ Ef{Ti)Ai{Ti) +0(h) + 0(d 4 



(*|P- ^H' g (M),Cg(Af , jV')J HtJ^) 

= {#,(>(), {H g (Af),H g (Ar)}} U ?(Tl)A , (Tl) + O(ft) + 0(rf 4 ) 



(9.12) 



(*|p- [c fl (^,jS0,c fl (^'>^")J n Tl |*> 



= {{iJ 9 (^),iJ 9 (Ar)},{iJ,(^) ) ^(^ / )}} U ?(Ti)jA , (Ti) + 0(^) + 0(rf 4 ) 



More generally, it can be shown that 



/ • \ n ~ 1 












a 


Y nj 


$3, 









n Tl |*) 



(9.13) 



{$ n , .... {$ 3 , {$ 2 , $!}} ...} L ?(TlW(Tl) + O(ft) + 0(rf 4 ) 



where each pair ($j,$j) can be set to be either the pair (i-L g (Aii), H g (A4i)), or the pair 
(Cg(M-i,J\fi), {Hg(Mi), H g (J\fi)}). In this sense the four- dimensional diffeomorphism algebra 
for pure gravitational fields emerges in the semi-classical limit for up to the corrections. 

The quantum spacetime including the gravitational and matter sectors, is invariant 
under the actions of four-dimensional diffeomorphisms acting on both sectors. On the other 
hand, the generators identified in (9.12) are expected to generate the four-dimensional dif- 
feomorphisms of the gravitational sector relative to the matter sector, and act nontrivially 
on the quantum spacetime |\&). It is these relative transformations that give the physics we 
observe. 
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10 Emergent Dynamics 

The Gauss constraint, momentum constraint and (modified) Hamiltonian constraint opera- 
tors are simply zero operators in the physical Hilbert space H. In this sense, the fundamental 
diffeomorphism symmetry becomes silent after providing the structure of the physical Hilbert 
space. Just as in classical general relativity, the physics in EI should to be instead given by 
relations between the dynamical fields. In the model, the relation between gravitational 
field and the matter coordinate fields are expressed by the local observables 0(X,T). The 
values {(fy\0(X, T)|^)} with different (X, T) are expected to be related through certain 
quantum transformations. Moreover, semi-classical limits of such transformations have clear 
classical counterparts that have to be agreed with; specifically, the transformations from 
(V\6(X,T)\V) to (V\6(X,T + AT) I*) has to give the clock-time dynamics of 0(X,T) in 
general relativity. In the last chapter, we found that the operators % g {M) and C g (M.,J\f) re- 
produce the algebra of relative diffeomorphisms between the gravitational and matter fields, 
and are candidates for generators of such transformations in the quantum theory. In this 
chapter, we explicitly evaluate the relative transformations of the emergent gravitational 
fields with respect to the clock time. We will show that these transformations are indeed 
generated by T-L' g (Af) and C g {M.,N), and they recover the dynamics of general relativity in 
appropriate semi-classical limits. 

Due to their self-adjointness, H g (M) and C g (M.,Af) generate the following unitary op- 
erators in A: 

exp (i 6 -±H g (M)) =t^,(jW,ci), exp(i^C g (M,Af)) = U Cg (M,^,e 2 ) 
When €\ is small, we have the infinitesimal transformations of the emergent gravitational 
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fields by Un^M,^) as: 



5 ei [_ EfiTjds, 



a 



e 1 (^\P i -[J(e XAX ) I ,W g (AOfeW + 0(e 1 2 ) + 0(h) + 0(d A ) 



(10.1) 



S ei Vexp[ 




(Ti)(Tj)de b }% 



ex,AX 



eimK[h(e XtAX )f , U g {N)]il Tl |*> + 0( £l 2 ) + 0(h) + 0(rf 4 ) 



Similarly, the infinitesimal transformations by Uc g (M.,N ,e 2 ) with small €2 are given by: 



Equations (10.1) are of special interest in deriving the dynamics in terms of the clock time. 
These equations are expected to give the variations of the gravitational fields at the clock time 
Ti, when the fields are translated perpendicularly to clock time slices. The perpendicular 
direction of the time slices might or might not agree with the clock time coordinate direction 
at the time T\. When the two agree, the clock time derivative of the local gravitational 
observables should be described by (10.1) with a specific choice of Af, and possibly a term 
generated by the Gauss constraint. In this case the matter coordinates simply correspond 
to the gauge where the shift function V a is set to be zero. When the two directions differ, 




e 2 (^\P l -[J(e XAX ) I , C g (M, M')]tl T M + 0(e 2 2 ) + 0(h) + 0(d 4 ) 




(10.2) 



6 2 (y\FUh(e XAX )lC g (M,M')]fl Tl \V) + 0(e 2 2 ) + 0(h) + 0(d A ) 
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the equations of evolution should further include a term that is generated by the momentum 
constraint. In short, aside from the contribution (10.1), the presence of the terms generated 
by momentum and the Gauss constraints should be determined by the dynamics of the 
spatial coordinate and frame fields. 

To calculate the clock-time derivative of a gravitational observable O, one needs to choose 
an appropriate value for M in (10.1). Ignoring matter back reactions, the commutator 
i[d(v*(r)),H- g (N)] approximates i[6(v*(r)),H(N)], which corresponds to {0(x), H(N)} 
in the classical theory. Recall that when all fields are on-shell, H(N = 1) generates 
the diffeomorphisms given by unit normal flows orthogonal to spatial slices, transforming 
0(x) = 0(x, 0) into 0(x, t). Therefore, when all the constraints are satisfied, (O(x), H(N = 
1)} = d t O(x,t)\ t=0 is simply the speed of O with respect to the proper time carried by the 
flow. Also, in chapter 5 we have defined the speed of the clock i^o(x) = {0° '(x) , H (N = 
1)} = dt4> (x,t)\t=o (on shell). Combining the two different speeds, Leibniz rule instructs 
us to obtain the derivative of O with respect to 0° at the proper time t — 0, by using 
v~o(x){0(x), H(N = 1)} = {0(x),H(N = z^o 1 )}- Therefore, we expect the appropriate 
lapse function to be Af(X) = v~o(X,Ti) . 

To independently verify these statements about the dynamics, we now carry out a cal- 
culation directly based on the definition of the observables. Since there is a set of local 
observables corresponding to each clock time, the dynamics can be evaluated by comparing 
the expectation values of the observables at a sequence of different clock times. In partic- 
ular, the comparison between two adjacent clock times leads to the equation of motion for 
(0(X,T)) as 



d 

dT 



(O(X.T)) = lrn. {*|PO(X)li, ; - j+A:/ -|*) (v]/|PO(X)llv,|*) (10.3) 

rp AT^O AJ 



To evaluate this quantity, we start from the invariance of the local observables. Being 
physical, a local observable 0(X,T) = FO(X)TL T commutes with H(N p ) for arbitrary N p . 
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Setting p = p* and N p = NiX we have 



= <tt|- \H (N$) ,PO(X)n Tl j |*) 

= (*\fo(x) 1 - \h (n$) ,n Tl j |*) + (*|pi [£ (N$),d(xj\ n Tl |*> 



(10.4) 



The second term is familiar from the previous chapter, while the first term can be calculated 
using the notation defined in (3.7) and (6.4): 



*|POpO- h(n p %),u Ti |*) 



(M/|p6(x)^| S[r] )( S[r] |^ «),n Tl ] |*> 



[r] 



£<*|pdpr)| S[F] ><^ [%) «) ,n ri (f)J p*//|*> 



^(*|PO(X)| S[r] ) 



[r] 



x 



[%)«), ^o(^(f))] v i (<(f))n^(<(f))5(0°(<(f)) -top^/i*) 



^rl^E [%)«)^°(C(f))] ^- n^«( f ))5(0°«(f)) - TOP*//!*) 

m m Tl „ 



+0(fc) 

(*|PO(x)J]jV(x m )d int (x ro )n ri |*) + (^\Fd(x)J2^(x m )o 4> o(x r , 



d 



dT„ 



n { r n} l*) 



Tl 



+0(d) + O(ft) 



(10.5) 



where ft{r„} gives the deformed spatial slices with nonconstant clock field values around Ti, 
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and is denned by 



( S[ r]|n { T„ } = (Sr\sym < i>#o(v*(r))5((j> (v*(r)) — T n ) > F di ff (10.6) 



for any (s[p]| G .ft. Since IT{ Tn=Tl } = 11^, the coherence condition (7.1) implies 



FO(X) 



d 



dT„ 



n { T n }|^) 



Tl 



d 



dT„ 



FO(X, {T n })fl Tl |^) + 



d 



Tl 



<9T„ 



10.7) 



popr,Ti)n {Tn} |*> + c>(ft) 



Tl 



where 0(X,T) can be any local observables involving only gravitational and the clock's 
momentum field, and 0(X, {T n }) is a smooth function with 0(X, {T n = TJ) = 0(X,T x ) = 
{0(X, Ti)). As proposed near the begining of this chapter, the lapse function will now be 
set to be Af(X) = v7 (X, Ti). Applying the coherence conditions to (10.5), we derive 



(®\FO(X) 



h 



H NJ 



(Ti) 



n 



0(X, Tk)(*|P v^ m , Ti)$ int (X m )n Tl |*> + 0(X, Ti) (*|P ]T 



g 

st: 



n { T„ } |^) 



Tl 



+O(X,T 1 )^^ 1 (X m ,T 1 ) 



9 



<9T„ 



^(X m ,{T n })(*|M Tl |*) 



Tl 



9 



m 

0(X,T x ) 



o(x, {T n })(*|pn Tl |*) + 0(d) + o(h) 



Tl 



^i/^(X m ,Ti) ($ int (X m ,Ti) + 



9 



<9T„ 



^o(A" m , {T n }) 



Tl 



d 

dT 



0(X,T) 



Tl 



+o(d) + o(n) 



(10.8) 



where we have used the ADM clock condition (6.1). The first term above involves the 
dynamics of the clock, and the second term is the desired clock time derivative of the local 
observable. 
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Further, combining (10.8) and (10.4), we obtain 



o = m 



7i 

= 0(X,T 1 



HlN™ ) ,PO(X)n Tl 



I*} 



d 



<9T„ 



v<p°(X m , {T n }) 



Tl 



n ri |#) + o(d) + o(ft) 



dT 



0(X,T) 



Tl 



(10.9) 



Specifically, when 0(X,T) = FU. T , equation (10.9) gives the clock dynamics 



^o(X m , Ti) ($ int (X m , Ti) + 



d 



dT„ 



po(X m , {T n }) 



Tl 



+ 0(d) + 0(h) 



(10.10) 



Finally, (10.9) and (10.10) leads to the relative Heisenberg equation with respect to the clock 
time 



d_ 
dT 



h 



+0(h) + O(d) 



0(X,T) 

6(x),h-(n p 



UtM + 0(d) + 0(h) 



lo.ir 



-i 

(Ti) 



n Tl |*> + (*|P 



0(X),^(iV p i Tl) 



n Tl |*> 



The second term in (10.7) involving H m requires some attention. Referring to (5.1) and 
(5.3), we recall that a spatially localized operator is a product of gravitational operators and 
the matter operators that provide the coordinates and frames. Accordingly the commutator 
between the spatially localized operator and H m (N p *) can be sorted into two parts. The 
first part comes from commuting the gravitational operators with H m (N p *), and gives the 
contribution of matter back reaction, which will be denoted as e m . The second part comes 
from commuting the matter operators with H m (N p *), and gives the dynamics of the mat- 
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ter coordinates and frames. From the previous discussion, this part should contain terms 
generated by momentum and Gauss constraints. 

Recall that the emergent gravitational fields transform as (8.2) between two spatial matter 
coordinates and frames. Consequently, we only need to evaluate the dynamics in one specific 
spatial coordinate and frame. Here we pick our specific coordinates and frames to be fixed 
at each v n , which means 



The dynamics of these specific matter spatial coordinates and frames have zero contribution 
to (10.11), up to an error of order of h. In this setting, the clock time derivatives of the 





(10.12) 
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gravitational local observables can be evaluated using only 7-L and the matter back reactions: 



d_ 
dT 
d 

dT 



(J( e X,AX,T) T ) 



Ti 



E?(T)ds a 



Ti •> S x ,ax 
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J(e XAX )j,H- g (N p 



"(Tl) 



n ri |*> + e m + 0(h) + 0(d) 



<tt|P- J(e XAX )i,W g (^o~ 1 (T 1 )) U Tl \^} + e m + 0(h) + 0(d) 
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dT 
d 

dT 



(h(ex A x,T^) 



Ti 



Vexplf A J b (T)(Tj)de b Yj 

Ti J ex, ax 



= (tf |P 



ft 



IL Tl \*) + trn + 0(h) + 0(d) 



(tf|P- ^(ex,Ax)J,^ nT 1 |^) + e m + 0(n) + 0(rf)) 



(10.13) 



By setting ei = AT and Af = i^o _1 (Ti) in (10.1), we see that (10.1) agrees perfectly with 
(10.13). This means that lig indeed generates the diffeomorphisms perpendicular to the 
equal-clock-time spatial slices, when acting on |^). 

Next, we evaluate the semi-classical limit of (10.13). The commutators in the equations 



68 



can be carried out again using (6.13) to obtain 



d , 



E?(T)ds a 



Sx,, 



= (*|p- [^.ax)/,^^- 1 ^))] n Tl |^> + e m + + o(d) 
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(*\$f AX ( AT,), HT,))^) + WS; a AX ( AT,), m), aiTj)^) + e rn + 0(h) + 0(d) 



dT^ 



Vexp 



A J b (T)(rj)de b 



ex, ax 



(*|p$f • AX (J, fe)£n ri |*> + e m + 0(h) + o(d) 

(*|^' AX (im), MTi))f |¥> + e m + 0(/») + O(d) 



;io.i4) 



Note that the equation for E fields contains a term $j a which involves the operator a. 
Finally, using the coherence conditions (7.1), we derive the equations of motion for the 
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emergent gravitational fields 
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+e m + 0(h) + O(d) 
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ex, ax 
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Pexp[/ A^r)^)^]^^^-^)) 

'ex ax 



B(Ti),i4(Ti; 



where H g (N) is exactly the classical gravitational Hamiltonian constraint. Note that the 
term of contributes only to the order of d 3 , and is suppressed similarly to the previous 
case. Since (10.15) holds for any Sx,ax an d ex,Ax, it gives the equations of motion for the 
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emergent fields up to errors of 0(d) 



d 

dT Ti 



Ef(X,T) = {Ef(X,T), [H g (N) + M g (V) + G g (A)}} 



E(Tj ) , A(Ti ) ,Af=^o - 1 (Ti ) ,V =0,A=0 



+e m + 0(H) + 0(d) 



_d_ 

dT 




(X,T) = {A t 



(X,T),[H g (Af) + M g (V) + G g (A)]} 



B(ri),A(Ti),Ar=i. o 



1 (Ti),y=0,A=0 



+e m + O(^) + 0(d) 



(10.16) 



Referring to (2.4) we conclude that (10.16) is the gauge-fixed (with N = z/^o -1 , V = 
and A = 0) classical equations of motion in Hamiltonian form, up to the corrections. Once 
again the correction terms are given by matter back reaction, quantum fluctuations, and the 
discretization of space. Up to these errors, with the matter coordinates and frames satisfying 
(10.8), the equations of motion agree with classical general relativity. 

The dynamics of emergent gravitational fields using other matter coordinates and frames 
follows immediately, from applying the clock time dependent transformation (8.5) to (10.16). 
Since (8.5) agrees with the classical transformations up to the errors, it will simply re-express 
(10.16) into another gauge in which V a and A 7 are nonzero. 

11 Corrections to the Classical Limit 

Before finishing our exploration, let us look into the correction terms that appear in the 
constraint equations (9.6) and (9.11), the emergent algebra (9.13), and the equations of 
motion (10.16). For the purpose of recovering general relativity, we assume that the matter 
back reaction e m is small and focus on other corrections. 

The corrections of order h result from two different kinds of quantum effects. The first 
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is the regularization of the inverse triad factor (det EY l l 2 appearing in H g . In the standard 
operator H g (N) defined in (2.29)-(2.31), this factor is regularized by the commutator between 
the total volume and holonomy operators, and thus is always finite by construction. This is 
true even when the state being acted on has zero total spatial volume. On the other hand, the 
classical counterpart (det-E) -1 / 2 clearly diverges when the E fields vanishes. This contrast 
signifies the departure of the quantum Hamiltonian constraint from the classical Hamiltonian 
constraint when the scale of space becomes comparable to l p . The gravitational Hamiltonian 
constraint operator H g (N p ) for our model, defined in (3.3)-(3.5), preserves this feature and 
therefore also contains the inverse triad corrections at the Plank scale. The second source 
of corrections is the uncertainty principle. All the fields are quantum mechanical in the 
model, so there are corrections coming from quantum fluctuations in both gravitational 
and matter sectors. In our context of describing local gravitational fields, the uncertainty 
in the matter fields results in fuzziness of the coordinates and frames. As a result, the 
conditions on the coordinates and frames listed in chapter 6 all come with errors of order h. 
Further, in obtaining the semi-classical limit, the replacement of the local gravitational and 
clock momentum observable operators with their expectation values also introduces errors 
of order hl 2 p and order h. Summarizing all the corrections denoted by 0(h), we see that they 
are negligible when the coherent state is observed at a large scale, when the competing 
classical terms are much greater. 

The corrections of order d result from the discretization of space. Remarkably, these 
terms are of zeroth order of h and could still dominate in large scales when the quantum 
effects are insignificant. Recall that d is the upper bound of the spatial coordinate gap, and 
it is finite since there is only finitely many physical points {v*} in the discretized space. This 
is in contrast with the regularization of a theory in a continuous space, in which coordinate 
gaps are merely regularization parameters to be set infinitesimal. The effect of the finite 
value of d comes in two parts. First, in the correspondence (8.1), we see that the expectation 
values of the local observables determine the emergent fields only up to errors of order d. 
These errors account for the fact that one can only specify a smooth field up to corrections 
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of order d, by making finitely many measurements on coordinate sites separated by gaps of 
typical value d. Indeed, this seems nothing more than a fact of reality, that we can only 
estimate continuous fields using finite data points. However, instead of the insufficiency 
of discrete data in describing a continuous reality, the model attributes the errors to our 
intention of obtaining a continuous picture from the intrinsically discrete reality. Second, 
the use of holonomies instead of the A fields as configuration variables introduces corrections 
in all orders of \A\d. In the regularized classical Hamiltonian constraint, the A fields appear 
through the linear terms in holonomies divided by the coordinate lengths of their paths. The 
nonlinear terms in the holonomies contribute to errors, and vanish in the limit of the paths 
shrinking to zero length. In the quantum theory, it no longer makes sense to take such limits, 
since the knot space is topological. In our model, we introduces matter coordinates such that 
we can define emergent A fields from the holonomy observables, through the correspondence 
(8.1). Similar to the classical case, when written in terms of the emergent fields, nonlinear 
terms of holonomy observables contribute to corrections in the semi-classical limit of the 
model. Moreover, the coordinate lengths of the embedded paths in the coordinate space are 
given by the finite coordinate gap d that can not be taken to zero as a limit. Therefore, 
the discretization of space in our model also leads to finite values of holonomy corrections, 
which contain all powers of \A\d. The holonomy corrections are clearly suppressed by small 
d value when A fields are nonsingular and bounded. However, the corrections would become 
important when A fields become singular near, for example, an initial singularity or a black 
hole. Overall, we see that all the corrections denoted by 0(d) are negligible when the matter 
spatial coordinate gaps are small and the emergent gravitational fields are nonsingular. 

Lastly, all the corrections in the model are well-formulated, and can be explicitly calcu- 
lated. Moreover, the model is ready for the inclusion of emergent matter fields defined in a 
way analogous to (8.10), and also their full coupling with the emergent gravitational fields. 
The explicit expressions of the corrections will be in terms of the emergent fields and their 
quantum fluctuations. Therefore, the model serves as a promising testing ground for the 
emergence of loop quantum cosmology from loop quantum gravity. 
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12 Conclusions 

Let us briefly review the results of this thesis. The thesis is based on the current stage of loop 
quantum gravity, which results from a specific canonical quantization of general relativity, 
and possesses a rigorously constructed knot space that solves the Gauss and momentum 
constraints. The space encodes the Planck-scale quantum geometry of space in a background 
independent way that respect the diffeomorphism symmetry in general relativity. Due to its 
background independence and discrete structure, the theory's semi-classical limit has been 
difficult and unclear. However, the success of loop quantum cosmology provides evidence 
of the possible correct limits. In parallel with loop quantum cosmology, the model in this 
thesis proposes a concrete method to overcome the difficulties and obtain the semi-classical 
limit of the full theory. 

The model we have constructed shares the same knot space with loop quantum gravity 
coupled to matter fields, although it uses the modified graph-preserving Hamiltonian con- 
straint operator in constructing the physical Hilbert space. Strictly respecting background 
independence, the model utilizes its matter sector to provide the coordinates and frames 
to describe the local gravitational observables. In obtaining the physical Hilbert space, the 
model assumes the validity of the group averaging procedure described in chapter 4. The 
semi-classical limits given by the coherent state |^) agree with full general relativity, when 
the criteria discussed in chapter 11 are met so that the corrections are insignificant. Also, 
the corrections of order h and d have clear interpretations, as discussed in the preceding 
chapter. 

Finally, the model provides a set-up to explicitly calculate the correction terms for the 
dynamics of emergent gravitational fields. It is of great interest to see how these corrections 
behave near the initial singularity, and near the singularity of a black hole. In loop quantum 
cosmology, the quantum and holonomy corrections have been extensively studied. Among the 
remarkable results of these models is that the corrections replace the initial singularity with 
a well-behaved bouncing of the scale factor, preceded by a contracting universe [31j[30j and 
followed by a built-in slow-roll inflationary phase [33J [34J [35J. These predictions are testable 
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and give distinguishable signals in the spectrum of cosmic microwave background radiation 
[33j[34][35j. It is then important to ask whether loop quantum cosmology emerges from a 
certain symmetrical semi-classical limit of loop quantum gravity, and if so, what additional 
details the full theory would provide regarding the predicted signals. Explicitly evaluating 
the correction terms in our model could provide answers to these questions. Hopefully, it 
would show that the observable predictions of loop quantum cosmology are indeed tests of 
loop quantum gravity, as a fundamental theory of gravity. (A concrete, affirmative result for 
emergence of the bouncing has been obtained and will appear in |36j.) 
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